NASA TN D-5906 


i 


* 


NASA TECHNICAL NOTE 



NASA TN D-5906 

a. t 


OJ 

XIRTLASD AFB, N M; uj 

Ll) 




PERTRAN - 

A TRANSPORT-PERTURBATION PROGRAM 


by John L. Anderson 

Lewis Research Center 
Cleveland, Ohio 44135 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


WASHINGTON, D. C. 


JULY 1970 


TECH LIBRARY KAFB, NM 




TECH LIBRARY KAFB, NM 


2733 


□ 13 


1. Report No. 

NASA TN D-5906 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

PERTRAN - A TRANSPORT- PERTURBATION 
PROGRAM 

5. Report Date 

July 1970 

6. Performing Organization Code 

7. Author( s) 

John L. Anderson 


8. Performing Organization Report No. 

E-5359 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 44135 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D. C. 20546 

10. Work Unit No. 

120-27 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

PERTRAN is written to be used primarily with the two-dimensional transport program 
TDSN. The equations of first-order transport perturbation theory are cast into three 
discrete forms consistent with the P-1, transport- corrected P-0, and diffusion approx- 
imations. The FORTRAN IV program calculates the neutron lifetime and effective 
delayed neutron fraction and the reactivity contributions of various reactions for each 
of the three approximations. Sample problems compare the three perturbation approx- 
imations with two-dimensional transport spatial calculations. Input instructions and 
listings of the program and a sample problem are provided. 


17. KeyWords (Suggested by Author(s)) 

18. Di stribution Statement 


Perturbation 


Unclassified 

- unlimited 


Computer program 





Nuclear reactor 





Reactivity 





19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

65 

$3.00 


*For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 


I 


■ ii in 






PERTRAN - A TRANSPORT-PERTURBATION PROGRAM 
by John L. Anderson 
Lewis Research Center 

SUMMARY 

PERTRAN is written to be used primarily with the two-dimensional transport pro- 
gram TDSN. The principal steps in the derivation of the equations of first-order trans- 
port perturbation theory from the Boltzmann transport equation are given. These equa- 
tions are cast into three discrete forms consistent with the P-1, transport-corrected 
P-0, and diffusion approximations to transport theory. The equations have been incor- 
porated into a FORTRAN IV program which will calculate the neutron lifetime, the effec- 
tive delayed neutron fraction, and the reactivity contributions of various reactions for 
the three approximations. 

One advantage of having three approximations is that by comparing the various 
approximations and their dependence on parameters such as mesh spacing one may 
choose an approximation for which the associated transport calculations require the 
least computer storage and computational time. Sample problems compare the three 
perturbation approximations with two-dimensional transport spatial calculations. The 
input and output features of the program are described, and listings of the program and 
a sample problem are provided. 


INTRODUCTION 

The design of a nuclear reactor requires many computer calculations to determine 
the nuclear characteristics. These calculations are usually made by one of the several 
multidimensional, multigroup transport and diffusion theory programs that are available. 
However, these programs often require large amounts of computational time. Further- 
more, the number of calculations needed is large because many design alternatives must 
be compared. Also, the effect of manufacturing tolerances and experimental uncertainty 
in the input cross sections must be determined. 

However, these design problems can be solved with greatly reduced computational 
time through perturbation theory. Somewhat fewer transport calculations, each of which 


is an iterative time-consuming process, are needed when the perturbation method is 
used. In using the method one first chooses a certain assembly as an unperturbed base. 
Then by using the transport theory solutions for this single unperturbed assembly one 
can determine the reactivity effect of small perturbations in the assembly. Perturba- 
tion theory for fast-neutron critical systems is described in reference 1. 

This report describes the computer program PERTRAN, which uses first-order 
transport perturbation theory to compute reactivity, neutron life-time, and effective 
delayed neutron fractions. PERTRAN is written to be used primarily with the TDSN 
transport program (ref. 3). Many perturbation programs have been written for both 
diffusion and transport theory (e. g. , ref. 2). PERTRAN differs from other transport 
perturbation programs in that it provides three approximations to the perturbation cal- 
culations. The best approximation uses P-1 cross sections and unperturbed real and 
adjoint fluxes and currents from the transport theory spatial calculation. The next best 
approximation is provided by transport- corrected P-0 cross sections and fluxes and 
currents. The diffusion approximation, which is least accurate, uses the transport- 
corrected cross sections and the fluxes from a transport calculation. 

These approximations allow one to choose the accuracy of the transport calculations 
to be commensurate with the importance of the perturbation. 


SYMBOLS 


B buckling factor 

E neutron energy 

H buckling dimension 

AAA 

i, j,k unit vectors in the x-, y-, and z-directions for rectangular coordinates (fig. 1) 
J(r,E) neutron current (directional) at position r with energy E 
k neutron multiplication factor; 1/k is the eigenvalue of Boltzmann equation 

1 neutron lifetime 

N q atomic density, atoms/(b)(cm) 

NG number of discrete energy groups 

NIJ number of discrete volume elements (V) 

P^ (p) Legendre polynomial 

P™( j u) associated Legendre polynomial 
P-0 zero-order ( l = 0) cross sections 
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* 

P-0 

P-1 

AQ 

r 

V 

x(E) 


y 

s 

6 

M 

u(E) 

$(r, E,ft) 


<?(r,E) 

* 

ft 


transport- corrected P-0 cross sections 
first-order (l = 1) cross sections 
incremental change in quantity Q 

position variable representing three-dimensional coordinates of neutron 
volume element, cm 

fission spectrum; probability that neutron released through fission will have 
a particular energy E: f x(E)dE = 1 

delayed neutron fraction of j n delayed group 

extrapolation distance constant (0. 71045608) 

macroscopic cross section, cm”* 

a, 

angle between k and ft (fig. 1) 
cosine of 8 

angle between ft and ft* 

average number of neutrons with energy E released per fission 

directional neutron flux; number of neutrons of energy E at position r 
flowing through a unit solid angle and unit area in direction ft 

scalar (nondirectional) neutron flux of energy E at position r 

A A ^ 

angle between i and projection of ft in the plane perpendicular to k 
direction of neutron flow at position r (fig. 1) 


Subscripts: 


g index of energy groups 

i index of spatial position 

l order of Legendre polynomial 

Superscripts: 

p perturbed quantity, s Q + AQ 

0, 1 order of Legendre polynomial for cross-section expansion 

t adjoint quantity 
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Cross-section definitions: 

(Macroscopic definitions are presented in discrete form for energy group g; the 
equivalent continuous form is shown for the second definition only.) 




2 

2 


N2N 

g-g* 


2 


( 0 ) 

g~g' 


2 


(1) 

g-g' 


2 


( 0 ) 

gg 


2 


Ik 


g 


2 


r 


g 


2 


S g~g' 


2 


2 


tr 


g 


diffusion coefficient, 



absorption (includes capture and fission), 2 (E) 

cL 


(Fick's law) 


fission 


n — 2n scattering from group g to group g r 

P-0 scattering from group g to group g' (includes elastic, inelastic, and 
twice the n — 2n scattering) 

P-1 scattering from group g to group g' 


within group scattering (used to provide neutron balance), 


40) s 

J gg 


S t or S tr 

L g 


2 - 2 
gj a g r g 


2 t for P-1 


2 tr for P-0 


transverse leakage (buckling- loss) cross section 
removal or outscatter, 2 r = ^ ( S g-g T " S g-g’) 


total scattering, 2 t 


"g-g' 


•g 

g#0 

= 2 ° , + 2 ? 


J g~g' g-g’ 


total, 2 t -2 +2 + 2^ 

g g g ^ 


transport, 2^ = 2 t - 

g g 


E 

g 


g 


TRANSPORT PERTURBATION EQUATIONS 

The time independent Boltzmann transport equation may be written 
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[£2 • V + 2 t (r,E)]#(r,E,£2) =-i- ff dE' d£2'[$(r, E', £2>(r, E’)S f (r, E’)x(E)] 

'ff dE* d£2* $(r, Ejfi’JSgOp, e' -* E,£2' — £2) (1) 


, a 

The continuous variables r, E, and £2 represent the dependence on position, energy, 

As 

and direction; <&(r, E, £2) is the real flux. 

Two other equations, nearly identical to equation (1), are needed to develop the 
perturbation equations. One equation provides the adjoint flux <3?^ (r, E, £2): 


[-£2- V + S t (r,E)]# 1 'Cr,E,£2) ff dE’ dfi’^ff, E',£2>(r, E)S f (r, E)x(E’)] 

47rk 1 ' 


ff dE’ d£2' (r, E ', £2')S g (r, E — E*,£2 — £2*) 


( 2 ) 


where = k. The other equation provides the perturbed flux <f> p (r, E, £2): 

[£2 • V + 2 p (r,E)]$ p ff,E,£2) = — ff dE' d£2* $ p (r , E ’ , £2' ) 

4nk p 


[^(?,E')2 f (?,E')] P x(E) + ff dE' d£2’ $ p (r,E’,£2')S p (r, E' - E,£2 f - 


£ 2 ) 


( 3 ) 


where each perturbed quantity is equal to the unperturbed quantity Q plus its per- 
turbation increment AQ (not necessarily small). 


Derivation 

The following derivation of the transport perturbation equations is taken from refer- 
ence 2. Some intermediate steps omitted in this report may be found in that reference. 

The transport perturbation equations may be obtained by multiplying equation (2) by 
4> p (r, E,£2) and equation (3) by 4>^ (r, E,£2), integrating the two equations over all space, 
energy, and direction, and then subtracting the resulting equations. The exact equation 
for an eigenvalue increment resulting from a perturbation is then 
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T + F + S 


( 4 ) 


±_ _ 1 

k P k 



P 


where 


-/// dr dE dQ, AZ t (r,E)4>^(r,E,&)$ p (r,E ,S2) 

F = — — Iff dr dE dE’ A|>(r, E')S,(r,E’)] x(E)$f (r, E)$ p (r, E’) 
4nk JJJ 1 

-///// dr dE dE’ dS7 d£2 r A[2 s (r,E f - E,$V - fi)]$^(r, E, fi)$ p (r, E\ ft') 


(5) 

( 6 ) 
(7) 


’■*/// dr dE dE>(r, E’)S f (r, E)] p x(E)<pt(r, E)<p p (r, E’) (8) 

Equation (4) determines the eigenvalue change resulting from a perturbation. The 
reactivity, however it is defined, can be subsequently obtained from the eigenvalue 
change. 

A first-order approximation to perturbation theory is obtained by replacing the per- 
turbed flux <3> p by the unperturbed flux 4>. This now restricts the perturbation to one 
which causes only a negligible change in the unperturbed flux. 

The unperturbed angular flux is then expanded in terms of spherical harmonics 
before performing the solid angle integrations in equations (5) and (7): 



l 

+ y ] (<P™ sin mi// + cos mi//)P™(jU) 

m-1 


(9) 


where 


ft = i sin 6 cos xp + j sin 0 sin xp + k cos 6 


and 


Id = cos 6 

(see fig. 1). If terms are retained only through the P-1 approximation, 

4?(r, E,ft) — — \cp (r, E) + 3ft * J(r, E)] (10) 

An 
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z(k> 



Figure 1. - Coordinate systems. 


where cp(r, E) is the neutron scalar flux and J (r , E) is the neutron current. A similar 
expansion for the adjoint yields 


4>^(r, E,0) (F, E) - 3f2 • (r,E)] (11) 

4n 

Furthermore, the scattering is assumed to be dependent only on the angle between direc- 
tions Cl and f? , and the incremental scattering cross section is expanded in terms of 
Legendre polynomials. The expansion is truncated after the P-1 term to yield 


AS (r, E* - E,n' - n) ss-i- M^E’ - E) + 3 u AS^V,E' - E )1 ( 12 ) 

s 47 r l ° J 

Substituting equations (10) and (11) in equation (5) and equations (10) to (12) in equa- 
tion (7) yields after integrating over the solid angle 


T = 



dr dE — AS t (r, E)[<?(r, E)^ (r, E) - 3J(r, E) • J(r, E)] 
477 L 


(13) 


and 


-/// 


dr dE dE’ — 
47T 


[AS^(r,E f - E)^,E>*(r,E) 


- SAS^^E’ - E)J(r,E*) • Jff,E)] (14) 
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Discrete Form 


The components of equation (4), which are equations (6), (8), (13), and (14), may be 
written in discrete form by replacing the continuous variables r and E by the discrete 
indices i and g. Hence, the equations may now be written as summations over volume 
increments and energy groups: 


NIJ NG 


T = > \ — AS. (P <p* - 3 J • jt V- 

ZZ zZ 4n V L 6 6 g gJi V i 


(15) 


i=l g=l 


NIJ NG NG 


F = 


z zztk a Ky . vgi^t 


(16) 


i-1 g=l g'=l 


NIJ NG NG 


S = 


- Z Z - 3as "'-gv • v » 


(17) 


i=l g=l g'=l 


NIJ NG NG 
V' 


p = ZZZiK']fwii v i 


( 18 ) 


i=l g=l g'=l 

This last equation is the product of the perturbed real power and the equivalent quantity 
from the adjoint calculation (an adjoint power). 


Approximations 

Equations (15) to (18) constitute the P-1 approximation to first-order transport per- 
turbation theory. Further approximations to equations (15) and (17) can provide a 
transport- corrected P-0 and a diffusion approximation. But before proceeding to these 
approximations, it is useful to isolate the reactivity contributions due to various reac- 
tions and regroup the P-1 equations (eqs. (15) to (17)) as reactivity sources or losses. 
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Transport - P-1 


The incremental total cross section may be written 


NG 


AS t„ - - E i2 g-g' - * As r g * * 41 ’ 


g g 

The reactivity contributions are as follows: 
Source: 

NIJ NG NG 


(Fission) F - - 


£££MhHv 4 v * 

i=l g=l g'=l 1 


(19) 


(20) 


NIJ NG NG 

(Scattering) S y y V' 

4tt / J J 


i=l g=l g T -l L 


AS^ ,(p . (p t, - 3AS^ t 

g-g’^gi^g' g 


^ t J • jt] V. 
-g g gJi 1 


( 21 ) 


Loss: 


NIJ NG 

(Absorption) A = j- g £ ^^4 ' 33 g ' 4 V i 


(22) 


NIJ NG NG r , 

(Removal) R = - £ £ ^ (AS g 2g'>. [v g 4 ~ V i 


(23) 


By collecting the current weighted terms in equations (20) to (23) we can estimate the 
contribution of the nontrans verse leakage out of the system: 


NIJ NG 


(Leakage) L = - — y y 
4 tt M ^ 


i=l g=l 


NG 

AS, J • jT - V AS,! 1 / „J„, • J 

V g g Lj g’-g g' g 
g g'=l 


V, 


(24) 


Note that the AS^. term includes the within group scattering; it does not cancel when 
current weighted as it does when flux weighted (eqs. (21) and (23)). 
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There is another reactivity loss mechanism which has not been accounted for - the 
transverse leakage or buckling loss. The cross section for the buckling loss, as used in 
reference 3, is 


S 



B 2 £. 


tr 


g 


( 2 tr « ♦ 2r) : 


(25) 


The buckling factor B is n/V3 for plane boundaries; for cylindrical boundaries, B is 
2(2. 405/ V3) and the buckling dimension H is the diameter. Within PERTRAN the in- 
crement AS., may be obtained by changing the transport cross section 2. or the 

^g g 

buckling dimension H. 

The contribution to the eigenvalue increment is 

NIJ NG 

<26) 

1=1 g=i 1 

which is identical in form to the equation of the total cross section (eq. (15)). If trans- 
verse leakage does occur, its contribution to the eigenvalue increment must be included 
in equation (4), which becomes 



F + S + A + R + TL 
P 


(27) 


Transport - Tran sport -Corrected P-0 (P-0*) 

If equations (15) and (17) are combined, 
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, NIJNG 
T + s=±y y 

4 * k h 


AS t g - 3j g * Jg ) - Z 

? ,= i 


NG 


( 0 ) 

g'-.gTg’Tg 


g 


NG _ 

+ 3 y AS^V _J_, • j t 
1-1 g -g g g 
g'=l 


i v ‘ 


NIJ NG 



K • ^ 


AS. 


NG m - ' 

V Asty J > 
L > g -g g 
'=1 


g 


'g 


+ (AS + 

a g 


aS r ) ivl) - Z iS g°-g v4 

g g'=l 


NIJ NG 

bEE 


in ■ 


i=l g=l 


tr g ^ 


g 

NG 


-3AS + „ J„ • J* + (aS^ + AS r ^<p a <p* 


i V i 


r g/ g g 


- T AS S'Pr'Vl 
g -g g g 

g’=l 


v. 


(28) 


where the incremental transport cross section is defined as 

N ° (1) 


iS tr = A \ - 

g g J 


y as j , 

g -g g 


(29) 


g 


Note that currents are needed from the spatial calculation in order to weight the trans- 
port contribution to the reactivity. The TDSN program calculates currents from the 
angular fluxes and, hence, can provide currents when using only P-0 cross sections. 
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Diffusion - P-0* and Fick's Law 


It requires a considerable amount of computer storage to provide both fluxes and 
currents for a spatial calculation. One further approximation is to use the diffusion 
theory definition of current, that is, Fick's law: 

j g = -D g *„ g (30) 

where is the diffusion coefficient for the group g. With this approach the currents 
can be calculated within the perturbation program from the gradients of the fluxes. 

The transport cross section term in equation (28) becomes 


3A2 j g ■ jt , 3D l AZ 

& e 


Vcp ■ V<? 


(31) 


If the fluxes <p come from a transport calculation, then equation (31) can be written in 

& 

terms of the transport cross section: 


3A2, ? ■ j£ = 

tr g g g 


1 ASt r (r 
1 g 


( 2 trj 


2 
S' J 


v«, g • V«,T 


(32) 


Equation (31) can also be written in terms of the diffusion coefficient: 

/ - aD e V 


3 A 2, 3 ■ jt = 

g g 


AD_ 

1 + — s 


X ■ H 


(33) 


The current terms in the transverse leakage contribution (eq. (26)) are also calculated 
using Fick's law. 

The program uses curve fitting techniques to determine the flux shape and subse- 
quently the gradients of the fluxes. The two particular techniques are presented in 
appendix A as they apply to fitting the flux profile to a second-degree polynomial. 


Prompt Neutron Lifetime 

The lifetime (l) of prompt neutrons (from ref. 2) is given by 
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where 


l = fff ^ dE E ’ <Fi E, &) 


(34) 


N = JJJ dr dE dE* x(E)^(r, E02 f (r, E')<p(r, E^ (?, E) (35) 

and (1/v) is the spectrum averaged inverse of the neutron speed. In discrete form these 
equations become 


NIJ NG 



Equation (37) for N is just the product of the unperturbed real and adjoint power - 
analogous to equation (18). In fact, in the absence of a perturbation to the fission cross 
section, the equation for the lifetime (eq. (36)) is identical to the equation for the reac- 
tivity contribution of an absorption perturbation (eq. (22) divided by eq. (18)) with AS 
replaced by (l/v) g . s 

Hence, the absorption cross section from a pure 1/v absorber can be treated sim- 
ply as an absorption perturbation. The perturbation in this case must extend over the 
whole region that established the flux spectrum A normalization factor such as an atom 
density N Q can be used so that the cross sections supplied are N Q (l/v) and thus the 
lifetime obtained is N Q Z . 


Effective Delayed Neutron Fraction 


The ratio of the effective delayed neutron fraction p ^ to the true delayed neutron 
fraction /3 for a delayed neutron group j is given (ref. 4) by 


hm Jff dr dE dE' i/2 f (r, E)<p (r, E)/3.(E’)</?t(r, E*) 
* fff dr dE dE' i/Ej(r, E)<p (r, E)x(E')<P^(r, E') 


(38) 
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where /3. (E) is the delayed neutron spectrum for delayed group j (normalized to /3j) and 
the denominator is simply the N of equation (35). 

In multigroup notation for /3- ( e ffy 



where g^ and g^ are the inclusive high- and low-energy groups that bound the partic- 
ular delayed spectrum j3- . 

V 

DISCUSSION OF SAMPLE PROBLEMS 

In order to determine the accuracy of the various approximations within PERTRAN, 
two-dimensional spatial calculations in x-y geometry were performed with the TDSN pro- 


^ Perfect reflection boundary condition 


1 - 


10 - 


15 — 


20 - 


Fuel (unperturbed) 



- x(cm) 
Perturbed 

kre& a f U ei regions 



-Molybdenum 

reflector 


'"‘“No return current 
boundary condition 


y(cm) 

Figure 2. - Sample problem geometry with boundary conditions and material re- 
gions (height = 45 cm). Perturbed regions for both sample problems are shown. 
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gram. These spatial calculations provided both the change in eigenvalue A(l/k) result- 
ing from a perturbation and the fluxes and currents which PERTRAN uses to determine 
this eigenvalue change. 

A quadrant of the fast spectrum reactor model is shown in figure 2. The fueled 
region is 30 by 30 centimeters in cross section; the annular molybdenum reflector raises 
the overall cross section to 40 by 40 centimeters. The reactor height is 45 centimeters. 
Four group cross sections from the GAM program (ref. 5) were used; the energy group 
structure is shown in table I. 


TABLE I. - ENERGY GROUP STRUCTURE 


Group 

Low energy boundary 

1 

a 0. 821 MeV 

2 

. 183 MeV 

3 

40. 87 keV 

4 

. 414 eV 


a Upper energy boundary is 
14. 9 MeV. 


Two sample problems were considered: (1) the perturbed region is small, at the 
center of the reactor, with a small transverse leakage rate; and (2) the perturbed region 
is larger, with a high transverse leakage, at the corner of the reactor. The pertur- 
bation in all cases was a 1-percent increase in material density in the perturbed region. 
The output listing of the center-perturbed sample problem is given in appendix B. 

For each of the sample problems the perturbed and unperturbed eigenvalues were 
calculated with TDSN using both P-1 and P-0 cross sections. The three approximations 
to the change in eigenvalue A(l/k) that PERTRAN provides are compared to the TDSN 
eigenvalue changes. 

In the absence of P-1 cross sections, TDSN will compute the current directly from 
the angular fluxes. Hence, the P-0 spatial calculation provided fluxes and currents for 
the corresponding approximation in PERTRAN as well as the fluxes for the diffusion 
approximation in PERTRAN. Furthermore, P-1 cross sections were used only for the 
perturbed material - not for the entire assembly. 


Center Perturbation 

For the smaller center region, the perturbation changed the eigenvalue 1/k by 
about 0. 005 percent (table II). All the PERTRAN approximations gave values of A (1/k) 


15 




TABLE n. - COMPARISON OF TRANSPORT AND VARIOUS PERTURBATION CALCULATIONS 


[Perturbation is 1-percent increase in material density; perturbed regions are corner and center (fig. 2); K e ^ f is nearly unity. ] 


Program 

Quantity 




Approximations 







Corner a 


Center (1-point) 13 


Center (3 -point) 0 



P-i d 

P-0 e 

Diffusion* 

P-1 

* 

P-0 

Diffusion* 

* 

P-0 

Diffusion* 

TDSN g 

k (regular) 

k (adjoint) 
k perturbed 

A(l/k) 

0. 999267 
. 999266 
. 999609 
-.342X10" 3 

0. 999508 
. 999507 
.999854 
-.346X10" 3 

0. 999508 
.999507 
.999854 
-.346X10 -3 

1.006531 
1.006531 
1.006586 
-.55X10" 4 

1.006541 

1.006541 

1.006595 

.54X10" 4 

1.006541 
1.006541 
1. 006595 
-.54X10" 4 

1.006262 
1. 006258 
1.006317 
-. 55xlO“ 4 

1.006262 
1.006258 
1. 006317 
-.55X10" 4 

PERTRAN 

A(l/k) 

Leakage 11 

-0.352X10 -3 

-.999xl0" 4 

-0.351X10" 3 

-.985X10 -4 

-0.276X10' 3 
-. 199X10“ 4 

-0.547X10" 4 
-. 102xl0" 7 

-0.548x10" 4 
-.102X10" 7 

-0. 551X10" 4 
-.391X10" 6 

-0. 550xl0" 4 
-. 140X10" 7 

-0.551X10" 4 

-.621X10" 7 


a Refers to perturbation of large corner region (see fig. 2). 

^Refers to perturbation of small center region (see fig. 2) using 1 mesh point in that region. 
c Same as footnote b but using 3 mesh points in that region. 

^P-l cross sections used in TDSN spatial calculations; corresponding approximation used in PERTRAN. 

e P-0 transport corrected cross sections used in TDSN; currents obtained and used in corresponding PERTRAN approximation. 

*Same TDSN calculation as in footnote e; TDSN currents not used in PERTRAN (KD--1 option used). 

^Convergence criterion used in TDSN was 10" 

^Leakage is in direction of calculation (not a buckling loss). 

that are within 2 percent of the TDSN value. The leakage (in the direction of the calcu- 
lation - in the xy-plane) was a negligible part of the total A(l/k) - about a thousand 
times smaller. 

5 (£ 

The P-1 and P-0 approximations provided nearly identical results for A(l/k) and 
the leakage. However, the leakage calculated in the diffusion approximation is a factor 
of 40 times greater than the P-1 and P-0 leakage. In this particular problem only one 
mesh interval was used in the 1-centimeter perturbed region. The same problem was 
also examined with three mesh intervals in the 1-centimeter region; the total number 
of mesh intervals remained the same (12 by 12). 

The leakage in the xy-plane was strongly affected by this mesh change. The P-0 
leakage increased in absolute value by about 40 percent; the diffusion leakage decreased 
to one-sixth of the value obtained when using only one mesh interval. These changes 
brought the diffusion leakage to a factor of 4 greater than the P-1 or P-0 leakage. The 
eigenvalue increments A(l/k) in the P-0 and diffusion approximations were essentially 
not affected by this mesh change. 

In order to further improve the leakage calculation (in the xy-plane) in the diffusion 
approximation, one would probably have to further increase the number of mesh inter- 
vals in the perturbed region, which in this case would then require more total mesh in- 
tervals. 


Corner Perturbation 


The corner perturbation produced an eigenvalue change of about 0. 035 percent. The 

* 

P-1 and P-0 approximations were the same, about 3 percent different from the TDSN 
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values of A(l/k). But the diffusion approximation to A (1/k) was only about 75 percent 

of the TDSN value. The P-1 and P-0 leakages differed by about 2 percent; the diffusion 

* 

leakage was about one -fifth of the P-1 or P-0 leakage. 

It appears that for the same accuracy the diffusion approximation will require more 
mesh intervals than the P-1 or P-0 approximations. However, if nontransverse leakage 
is an important part of the total eigenvalue change, then the diffusion approximation may 
not be adequate regardless of the number of mesh intervals. 

The PERTRAN input instructions and notes on the input and output features are pre- 
sented in appendix C. Appendix D is a FORTRAN IV listing of the entire PERTRAN pro- 
gram. 


CONCLUDING REMARKS 

The most important characteristic of PERTRAN, compared to other perturbation 
programs, is that it offers three approximations to the perturbation calculation. The 
approximations, P-1, P-0 , and diffusion, are in order of decreasing accuracy, but the 
associated transport calculations also require decreasing amounts of computer storage 
and computational time. 

The most accurate (P-1) approximation requires the use of P-0 and P-1 cross sec- 

* 

tions and, hence, the calculation of currents as well as fluxes. The P-0 approximation 

5): 

in PERTRAN requires only a P-0 cross section set but it also requires the currents 

* 

from the spatial calculation. The diffusion approximation (least accurate) uses only P-0 
cross sections and fluxes. (The fluxes, however, may come from a P-0 spatial calcula- 
tion and so they do have that accuracy. ) The leakage is obtained from the flux gradient 
(using Fick's law). Since only fluxes from the spatial calculation are required, this 
approximation could treat problems with many more mesh points and groups than the 
approximations requiring currents. 

The relative accuracy of the approximations and the dependence of accuracy on num- 
ber and spacing of mesh intervals may be determined for a specific problem. Thus, an 
approximation which gives a certain accuracy for the least amount of computer storage 
and computational time for the spatial calculations may be used. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 28, 1970, 

120-27. 
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APPENDIX A 


CURVE-FITTING TECHNIQUES FOR FLUX PROFILE 


The first technique, the Vandermonde matrix method, which will fit a polynomial of 
degree n-1 through n given points, is used in PERTRAN to obtain the flux shapes and 
then the gradients. Three points or fluxes at three consecutive mesh intervals are used 
to fit a second-degree polynomial. The following derivation of the equations is con- 
ducted in generalized form in which y will represent the flux and x will represent the 
spatial variable. The polynomial is 

2 

y = a Q + ajx + a g x (Al) 


which is Xa = y in matrix form. The column matrices (a and y) are 


a 



\ 

/ 


y 



and the Vandermonde matrix is 


h 


X = 


1 


V 1 



(A2) 


(A3) 


(A4) 


If X is written as the product of a lower triangular matrix L and an upper triangular 
matrix U, then 


X = LU 


(A5) 


and the inverse of X is given by 



X -1 = IT 1 !/ 1 


(A6) 


■with 


U 


1 ~ X 1 x l x 2 

0 0 -x^-Xg 

0 0 1 


(A7) 


and 


( 


-1 


0 

1 


X 1 “ x 2 


x 2“ x l 


0 

0 


A 


V (Xj - X 2 )(x 1 - Xg) (x 2 - X 1 )(x 2 - Xg) (Xg - X^Xg - Xg )^ 


Now the matrix of the coefficients a j is given by 


a = X -1 y 


a = U _ 1 L _1 y 


Evaluating a gives 


-1 


yi 


yi , y£ 


A 


yi 


X 1 - x 2 x 2 X 1 


y2 


ys 


^1 - X 2 )(x 1 - Xg) (X 2 - Xj)(x 2 - Xg) (Xg - Xj)(Xg - Xg ) . 


(A8) 


(A9) 


(A10) 


Let the lower element of the Ly matrix be z. Then, 
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(All) 



x iyi 

(*1 - x 2 } 


x ^2 
( x 2 " x l) 


+ x l x 2 



U 



(xj + x 2 )(z) 


This is the matrix of the coefficients a^. 

The gradient of y at any location x- within the range of validity of the three-point 
fitted polynomial is given by the derivative of equation (Al): 


vy 


= dy 
X j ^ 


= yU = a i + 2a * x ; 


’X] 


2 J 


X i 


(A12) 


Substituting the coefficients a^ and a 2 from equation (All) gives 



1 2x- - (xj + x 2 ) 

« \7 . j,. 

l 2x- - (x x + x 2 ) 

X 3 

+ 

Xj - x 2 (Xj^ - x 2 )(x 1 - Xg) 

y l + 

x 2 " X 1 ( x 2 “ x l )(x 2 - x 3 } . 


2x- - (xj + x 2 ) 
( x 3 " x l )( x 3 “ x 2 } 


y 2 

• ys 


(A13) 


This is the desired flux gradient (Vcp) expression. Through Fick’s law the flux gradient 
determines the neutron current for a particular energy and location: 

J = -D Vcp (A14) 

However, because TDSN provides scalar fluxes for the midpoints of mesh intervals, 
the flux gradient at these midpoints does not necessarily provide a good representation 
of the current through that mesh interval. More information can be incorporated into 
the current calculation by using the net current through the parallel faces of the mesh 
interval to represent that interval. Thus, a difference in gradients taken at boundary 
points of mesh intervals is used as the net flux gradient Vcp : 

v y j+ l - v yj = 2a 2 (x j+l ' x j ) (A15) 
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Vy i+1 ■ Vy i r |( Ax 2 )yi " ( Ax l + Ax 2 )y 2 + (Ax l )y 3| (A15a) 

J J (Ax^ • AXgjiAXj + Ax 2 ) J 

with AXj = x 2 * x l anc * Ax 2 = x 3 " x 2 * 

Frequently, a calculational configuration has boundaries that are perfectly reflect- 
ing (i.e. , the cell condition). In such case there is no net current across the boundary 
and, therefore, the flux gradient is zero. If a flux profile fitted to three mesh points is 
extrapolated to the boundary, it will not necessarily have zero slope at the boundary. 
Therefore, if the cell condition exists, the PERTRAN program fits a second-degree 
polynomial to the slope (identically zero) at the boundary and the two closest flux points. 
Equation (Al) is the generalized polynomial for the flux profile. The flux gradient is 
then 


y' = a^ + 2agX (A16) 

The zero gradient restriction 

y b = 0 = a l + 2a 2 x b ( A17 ) 

at any interior or exterior boundary x^ determines a^ in terms of a 2 and x^. 
Therefore, using the flux (y^ and y 2 ) at two adjacent mesh points (xj and x 2 ) we 
obtain from equation (A16) 

a 2 = r - - Yl - -r (A18) 

( x 2 - x l> L (X 2 ' X b } + (X 1 ‘ x b>J 

Now to determine the net flux gradient across such a bounding mesh interval with a cell 
condition on one side, we need merely to evaluate equation (A16) for y' at the boundary 
between x^ and x 2 . This particular gradient is given by 

y* =2a 2 (x-x b ) (A19) 

and is used as V<p in equation (A14). 
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APPENDIX B 


lllllllllllllll llllll III 


SAMPLE PROBLEM OUTPUT 


This sample of PERTRAN output is the corner perturbation problem (P-1 approxi 
mation) discussed earlier. 


SAMPLE PRJBLtM CORNER P-1 APPROXIMATION 


9 MAT 

NU 

NF AST 

NDOWN 

1 

A 

A 

3 

KCTYP 

KIN 

KPAP 

KP 

0 

1 

0 

1 

N2JNE1 

NZ0NE2 

LI STFX 

NMID 

3 

3 

1 

3 

LAST 

LAST1 

LAST2 

XKEFF 

A 136 

A136 

AO 18 

0.99926700 


l» J ) 

O.AOOOOOOE 

01 

O.AOOOOOOE 

01 

O.AOOOOOOE 

01 

3 . A003000E 

01 

0.2CO0000E 

01 

0. 2 0000 OOl 

01 

U- 5 000 00 J E 

01 

0 . 5003003 E 

01 

C.ACOOOOOt 

01 

O.AOOOOOOE 

01 

O.AOOOOOOE 

01 

O.AOOOOOOE 

01 

C. 2000000E 

OL 

0.2 OOOOOOE 

01 

0. 5 OOOOOOE 

01 

3.5 OOOOOOt 

01 

o.AOoooooe 

01 

0. A OOOOOOE 

01 

O.AOOOOOOE 

01 

O.AOOOOOOE 

01 

C. 2000000E 

01 

0.2 OOOOOOE 

01 

0.5000000b 

01 

0.5 OOOOOOE 

01 

O.AOOOOOOE 

Ol 

O.AOOOOOOE 

Ul 

O.AOOOOOOE 

01 

O.AOODOOOE 

01 

0.2000000E 

01 

0.2 OOOOOOt 

01 

0.5 OOOOOOE 

01 

0.5003000E 

01 

O.AOOOOOOE 

01 

0. AOOOOOOt 

01 

O.AOOOOOOE 

01 

O.AOOOOOOE 

01 

0 . 2 00 U 0 OOt 

01 

0.2 OOOOOOE 

01 

0. 50000 JOE 

Ol 

0. 5 OOOOOOE 

01 

0. 20000 00E 

01 

0.2 OOOOOOt 

01 

0.2 OOOOOOE 

01 

0.2 003000 E 

01 

0. 1 OOOOOOE 

01 

0.1 OOOOOOE 

OL 

0. 2 5 000 00 E 

01 

3 . 2 503 003 £ 

01 

C.2C0U000t 

01 

0.2 OOOOOOt 

01 

0.2 OOOOOOE 

Ol 

0 . 2 003003 E 

01 

0. 1Q00O0OE 

01 

0.1C00000E 

01 

0.2 5 00000 E 

01 

3 . 2503000E 

Ol 

0. 2000000k 

01 

0.2 COOOOOE 

01 

0. 2 OOOOOOE 

01 

3.20000006 

OL 

c . iooouooe 

01 

0. 1 OOOOOOE 

01 

0 .25000 00 E 

01 

0 . 2503000 E 

Ol 

C. 2 COOOOut 

01 

0.2 OOOOOOE 

Ol 

0.2000000E 

01 

0.20000006 

01 

0. 1000C00E 

01 

0. 1 OOOOOOE 

OL 

0.2500000b 

01 

0.2 500000E 

01 

C.2C0OO OOt 

01 

0.2 OOOOOOE 

Ol 

0 . 2 0 00 0 00 E 

01 

3.2 000000 £ 

01 

C. lOuOUOOt 

01 

0. 1 OOOOOOE 

01 

0.2500000E 

01 

0 . 2503 000 E 

01 

C.5000000E 

01 

0. 5 OOOOOOt 

01 

0.5 OOOOOOE 

Ol 

0 .5 OoDOOO 6 

01 

C.250000UE 

01 

0. 2 500000E 

01 

0 .6 2 503 00 E 

Ol 

0 .6250000 £ 

01 

0.5G0OU00E 

01 

0. 5 OOOOOOE 

01 

0.5000000 E 

Ol 

0 . 5000 000 E 

01 

0.25000 OOt 

01 

0.2 500000E 

01 

0. 6250000E 

01 

3.6 250000E 

Ol 


NUP 

0 

KS 

1 

KAPROX 

2 

NDSP 

-0 

O.AOOOOOOE 01 

NR 

12 

KAR 

1 

KF LUX 
1 

1GHI 

-0 

0 . 20000 006 01 

NZ 

12 

KO 

1 

KBEFF 

•3 

I3L0 

-3 

0. 2 OOOOOOE OL 

3.23333336 

31 

O.AOOOOOOE 

01 

0.2 OOOOOOE 

01 

0. 20000 OOE 

01 

0.2333303E 

01 

O.AOOOOOOE 

Ol 

0.2000000E 

01 

0.2 00 00 OOE 

01 

0.20300006 

01 

O.AOOOOOOE 

01 

0 • 2000000E 

Ol 

0. 2 0000 OOE 

01 

0.2303000E 

01 

O.AOOOOOOE 

01 

0.2000000E 

01 

0. 2000000E 

31 

0.2303000E 

01 

0.2000000E 

01 

0. 1000000E 

01 

0.10000006 

01 

0. 1003003F 

01 

0.20030 006 

01 

0. 10000006 

Ol 

0.1000003= 

OL 

0. 1 033303 F 

01 

0.20000006 

01 

0. 1000 0 COE 

01 

0. IOOOOOOE 

Ol 

0. 13333036 

31 

0. 2000000b 

01 

0. 1000000E 

01 

0. 1000330E 

31 

0. 1333.303F 

01 

0.2000000E 

01 

0 . IOOOOOOE 

01 

0.1 OOOOOOE 

01 

3. 133000 3 F 

01 

0.50000006 

Ol 

0.2500000b 

01 

0.2500000= 

01 

0. 25333006 

01 

0.5000 000b 

01 

0.2500000E 

01 

0. 2 500300E 

01 

0. 2533003E 

Ol 


0. 

O.UOOOOOE Oil 


U. 20Q0U00E 01 
G.lAOOUOOt 02 


O.AOOOOOOE 01 
0. 1 5 00 0 00 E 02 


0 .6000000 E 01 
3 . 1 753000 E 02 


0. 8000000E 01 
0.2000000E 02 


0.1 000000E 02 


0. 11 00000c 02 


3. 1 203000F 02 


Z( J) 

0. 

C.1300000E 02 


0.2000000E 01 
O.IAOOOOOE 02 


O.AOOOOOOE 01 
0.1 5 000 00E 02 


0 .60000 00 E 01 
3.1 750300E 02 


0. 8000000E 01 
0.2000000E 02 


0.1000000E 02 


O.linoOQOE 02 


0.1203000E 0? 


DfcLRl I) 

0. 2000000E 01 

c.ioouoooe ol 


0.2 0000 OOt OL 
0. 1 75 OOOOt 01 


0.20000001: 01 

0. 2 5000 00 E Ol 


0 . 2000000 E 01 


0. 1500000E 01 


0.1 OOOOOOE 01 


0.1000000E 01 


O.L003900F 01 


UtLZI J ) 

C. 2000 0 OOt OL 
C.1C00000E Ui 


0.2C00000E 01 
O.M500006 01 


0.2 OOOOOOE 01 
0.2 5 000 00 £ 01 


O.20OUGOOE 01 


0 . L503000E 01 


0.1000000F 01 


0. 1 0000 OOF 01 


0. 1003000E Ol 


Chit I G I 

0 . 75 3 S2 C Ct CO 0. 21 Jb 700E 00 0.2963000E-0 1 0 . 28800006-02 


DELTAS COMPUTED INTERNALLY FROM UlFFERENct BETWEEN PERTURBED AND UNPERTJRBFO CROSS SECTIONS 
TOSN FJKMAT cRJSS SECTIONS. TRANSPURT APPROXIMATION. WITH P-0 AND P-l CS . 


C.l 7 39SS5E- 03 
0.195999NE- 03 
0.29599 9S£- U3 
0.53 19999t- 03 


0. A360U01E-03 
0. 3960002E-03 
0. 5 300003E-03 
0. 8AJ9999t**03 


MATERIAL 1 
0.127 1000E-02 
0. LA95000E-02 
0- 1968000E-02 
0.2A2900LE-02 


0 . 8 730 OOAE-O 3 
0.12A6999E-02 
0. L651000E-02 
0.1897000E-0? 


0. 

0.20199 99E- 03 
0 .AA999996-0A 
0.2000001E-0 A 


0. 

n. 

0. 2000001E -04 
0. 5 00000? f -35 


0. 

3. 

0. 

3. 


03D301F-35 


TOSN FORMAT. P-l CRuSb SECTIONS nl Trt FACTOR OF 3 INCLUDED 


MATERIAL 1 
0.8A10001E-03 
0.1 L35000E -02 
0.1 718000E-02 
0.232 90 01E-02 


0. 129A000E-02 
0. 1023000E-02 
0 . 7263 00 A E -03 
0 . 2023 00 1 E-03 


0. 

0. 

-0. 1999999E-05 
-0.2000001E-05 


DELTA UUTSCATTER CROSS SECTIONS 

0.2 2 5 11 OCE- 03 0. 506 6000E-OA 0.2002000E-QA -0. 


SUCKLING LOSS INFORMATION **** 

FIRST DIMENSION A5. 000000 SEC ONO DIMENSION -0. 

FIRST DIMENSION Pt K TURSED -0. SECOND DIMENSION PERTURBED -0. 

BUCKLING FACTOR 1.813799A 


UN.1G, ITYI FOR 1 TY OF SIG A, NU SIG F,ID* UR TR OR T OT J , SIG OJTSCATTER, NUP VALUES OF SIG UP, AND NDOWN VALUES OF SIG DOWN. 


0.173999 SE- 03 
0.195999 9E- 03 
0.2959999E- 03 
0.531 99 9 St- 03 


0. A360001E-03 
0. 3960002E-03 
0. 3300003b- 03 
0.8A39999E-03 


MATERIAL 1 
0 • 127 1000E-02 
0 . LA95000E-02 
0.1 968000 E— 02 
0.2A29001E-02 


0.2251 100E-03 
0.5066000E-0A 
0.2002000E-OA 
- 0 . 


0. 

0.2019999E-03 

0.AA99999E-0A 

0.2000001E-QA 


0. 

0. 

0. 2000001E-OA 
0. 5000002E-05 


0. 

3. 

0. 

0. A033331E-35 
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C(NflG'ITy) FUR ITY OF 0.0 , P-0 G-TO-G. TRANSPORT • P-1 G-TU-G, NUP VALUES OF StG UP, AND NO J«IN VALUES OF SIG D"HN. 


0. 

J. 8 730004E-03 

MATERIAL l 
0.841000 IE -0 3 

0. 1294000E-02 

0. 

0. 

0. 

0. 

0.1246999E-02 

0. 1 135300E.-02 

0. 1023000E-02 

0. 

0. 

7. 

0. 

0. 1651000E-02 

0. 1 71d000E-02 

0. 7260004E-03 

-0.1999 999E-05 

0. 

0. 

0. 

0. 1 89 7 000E-02 

0.2 32900 IE-02 

0 • 2 02 ) 00 IE - 0 3 

-0.2000001E-05 

0. 

0. 

DELTA BUCKLING LOSS CRUSS SECT1UNS 
FIRST DIRECTION 

-0.A62 101 2E- 04 -0.4918664E-04 

-0.461353 3E-D4 

-0.4080 136E-04 





the 

NON 

-ZfcRO 

MATERIALS 

1 ARE 

PERT URBE D 

0 


0 


l 








MAP 

MAC 1, 

J) 









1 

l 

1 

1 

1 

1 

1 

1 

1 

l 

2 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 

2 

1 

i 

1 

1 

l 

1 

1 

1 

1 

1 

2 

2 

1 

1 

l 

1 

1 

1 

1 

1 

1 

1 

2 

2 

l 

1 

1 

1 

1 

3 

3 

3 

3 

J 

2 

2 

1 

1 

1 

1 

1 

3 

3 

3 

3 

3 

2 

2 

1 

1 

1 

l 

1 

3 

3 

3 

3 

3 

2 

2 

1 

l 

1 

1 

1 

3 

3 

3 

3 

3 

2 

2 

1 

1 

1 

1 

1 

3 

3 

3 

3 

3 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 


MAP MA(I,J) kITH ONLY PERTURBED ENTRIES 
000000000000 
000000000000 


l I 1 1 1 


i 

.m.E- 01 *. ! - 


I i= 
J - 
13 = 


A!3jon;T FLUXED y::a(i,j) 
1 , 


THE 
I 3 = 


MRoT I ‘IRK'JTION RhKNTS XJ(I,J) 

1 , 


MRJT .■IRE'.TI'H Al’JoiliT -'Rh-MTS XJA(I,J) 
1 , 


THE 

13 = 


THE 
I i - 
J • 

;i'E 

j 

iHE 

J 


TE * n ! 

1 , 

1 , 

.‘E ' 

1 , 


i : IK** -:i >u<rNT.; 

;• : if- **■ <vj rrr v.* •% ( i , / 


l .-.KH :;'Y r 1 1 , j 

1 


.r i -M-h TY F A ( I , J ' 


■\ tCR * r- tv r hR “ . ’ A! 0. 


PE«TURB6Di-ll/KtFF UNPERTURBED)) 
l TRANSPORT APPROXIMATION 


PROUUC T IUN 
— 0.3102826E-0J 
-C.9C94B46E-C4 
-C. 130 75 62L-04 
-C.1202218E-G5 
GKUUP TOTALS 
-C.41550d9E-0J 

TRANSVERbt 

leakage i 

-0. U72d66fc-04 
-0. 1896996E-04 
-0.85212 7 BE-05 
-C. 10O4345E-C5 
GROUP TOTALS 
-0. 4022424E-04 


SCATTERI NG 

0 . 

-0.64783606-04 
-0.2598379E-04 
- U. 6 95 76 52 E- 05 

-0.9 7 72 50 At — 04 

TRANSVERSE 
LEAKAGE 2 

0 . 

0 . 

0 . 

0 . 

0 . 


{ U/KEFF 
P- 

ABSORPT tUN 
0 • A A 163161-04 
0. 75 59 Id 7t-0A 
0.5A67L 71E-0A 
0. 1 3 095 A 3 E-04 

0. 18752 22k -03 


REMOVAL 
0.5 7135 A8 E“0 A 
0. 1953820E-04 
0.3697 730E-05 
- 0 . 

O.BQ371AOE-OA 


G-TO-G***J-riGT 
-0.2A8S 1 l 3E-0A 
-0. 258154 OE -04 
-0. 128564 7E-0A 
-u. 2604 322E-05 

-0.6615832E-04 


L EAKAGfc 

-0. A 7421 ABE -04 
-0. 328345 7E -04 
-0. 16 1603 3F - 04 
-0. 3498 7 91 F - 05 

-0.9991522E-04 

GROUP 

-0. 245593 7F -03 
-0. 1053 884E-03 
-0.206771 3E-05 
0. 1 32689 3C - 05 

-0. 3517229E-03 


UNPERTURbfcO MULTIPLICATION FACTOR **** 0.99926700 

PtRTURBED MULTIPLICATION FACTOR **** 0.99961832 


APPENDIX C 


PROGRAM INFORMATION 
Input Instructions 

This section contains the input instructions and explanations of the input parameters. 
The symbol * after a card number means to use as much of the card or as many cards 
as necessary. 


Card Format Variable 


Description 


1 II, IX, 14A5 TITLE Title cards. Number in card column 1 signifies the 

last title card. 

2 7110 NMAT Number of materials for which perturbed cross sec- 

tions will be provided. 

NG Number of energy groups. 

NFAST Number of fast groups (to be used if KCTYP - 0 or 1 
and NUP > 0). 

NDOWN Maximum number of groups down- scattered to. 

NUP Maximum number of groups up-scattered to. 

NR Number of first direction mesh intervals (horizontal 

direction on map - left to right). 

NZ Number of second direction mesh intervals (vertical 

direction on map - top to bottom). 

3 7110 KCTYP Format for cross sections (see card 15). 

= 0 Perturbed then unperturbed cross sections in 
TDSN format. Increments are obtained within the 
program: AS = S^ - S. (This option must be used 
if a buckling loss is to be determined - card 16. ) 

= 1 Increments AS in TDSN format. 

KIN Format for geometric specifications. 

=1 Use mesh from TDSN binary dump (obtained from 
TDSN option KBCDUP = -2). See also cards 6 to 8. 
= 2 Use binary dump of Ar and Az instead of r 
and z. See cards 9 and 10. 

= 3 Not binary (see cards 11 to 13). 

KPAP = -1 Read in perturbed v - fission cross sections, and 
then the adjoint production (Pj ). 
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Card 


Format 


Variable 


Description 


4 7110 


KP 

KS 

KAR 

KD 


NZONE1 

N ZONE 2 

LISTFX 

NMID 

KAPROX 


KFLUX 


= 0 Read in real (P.) and adjoint (pt) production 


NIJ 


= 1 Read in PAP, which is ^ ^ P^P^V^ (cannot be 


used if KBEFF = 1). 

= 1 Contribution to A(l/k) from production incre- 
ments will be calculated. 

= 1 Contribution to A(l/k) from incremental scat- 
tering into a group will be calculated. 

= 1 Contributions to A(l/k) from transport, absorp- 
tion, scattering removal, and buckling loss incre- 
ments will be calculated. 

= 2 Lifetime will be calculated. 

= 1 Contributions to A(l/k) from diffusion coeffi- 
cient increment will be calculated. Use only with 
KAPROX = 0. If KCTYP * 0, AD will be calcu- 
lated internally as AD = A £l/ (32^ )J . Restricted 
to NR £ 3 and NZ =£2. 

= -1 Same as for +1 except that cell boundary condi- 
tion exists which will be specified on card 30. 

Number of material zones in first direction (corre- 
sponding to material map in TDSN but not restricted 
to that map). 

Number of material zones in second direction (corre- 
sponding to material map in TDSN but not restricted 
to that map). 

= 1 List flux input and production rate input as part 
of output. 

Number of materials in identification map (IDM) 
(hence, TDSN map can be used directly). 

Type of approximation. 

= 0 Diffusion theory. 

jje 

= 1 Transport theory - (P-0 ) transport corrected 
P-0 cross sections. 

= 2 Transport theory - P-1 cross sections. 

Format for fluxes and currents (see card 21 ). 

= 1 From TDSN binary dump. 

= 2 5(15, E10. 6). 

= 3 7E10. 6. 
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Card Format Variable Description 

KBEFF = 1 Perform /3-effective calculation; read in de- 

layed spectra information at cards 32 and 33. 

5 F10. 8, 3110 XKEFF Multiplication factor (k) from unperturbed calcu- 

lation. Must be included for all problems, even 
if KBEFF - 1. 

NDSP Number of delayed spectra to be read in. 

IGHI Highest energy group in which any of NDSP 

spectra contribute. 

IGLO Lowest energy group in which any of NDSP 

spectra contribute. 

If KIN = 1 and KBEFF = 0 read in cards 6 to 8. If KIN = 1 and KBEFF = 1 read in only 
card 6. 

jjc 

6 Binary V(IJ) Volumes from TDSN binary dump. 

j|( 

7 Binary R(I) Mesh boundaries in first direction (NZONE1 

values). 

jjc 

8 Z(I) Mesh boundaries in second direction if NZ > 2 

(NZONE2 values). 

If KIN = 2, read in cards 9 and 10. 

a|c 

9 Binary DELR(I) Mesh increments in first direction. 

j|e 

10 Binary DELZ(I) Mesh increments in second direction. 

If KIN = 3, read in cards 11 to 13. 

11 110 KGEO Geometry. 

= 1 Slab (x - y). 

= 2 Cylinder (r - z). 

= 3 Sphere (r). 

12* 5(15, E10.6) NM(I), RM(I) First direction mesh. NM is the number of mesh 

intervals to include between the preceding value 
of RM and the value of RM that immediately 
follows NM. If NM = 0, associated RM is 
ignored. If NM < 0, associated RM is the 
last value to be used. SUM(NM) = NR. 

13* 5(15, E10. 6) NM(I), RM(I) Second direction mesh if NZ > 1. Same as for 

card 12 except that SUM(NM) = NZ. 
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Card Format Variable Description 

If KBEFF = 1 read in only cards 22, 28, 29, 32, and 33. 

14* 7E10. 6 CHI(IG) Fission spectrum (NG values). 

5|c 

15 TDSN C Cross sections (see Input Notes, p. 29). 

7E10. 6 If KAPROX = 2 and KCTYP = 0, then the order 

of cross sections for each of NMAT materials 
is as follows: 

Perturbed P-0 c.s. 

Unperturbed P-0 c.s. 

Perturbed P-1 c.s. 

Unperturbed P-1 c. s. 

Perturbed removal c.s. 

Unperturbed removal c.s. 

If KAPROX < 2, do not include P-1 cross sec- 
tions . 

If KCTYP = 1, each set of perturbed and unper- 
turbed cross sections is replaced by one set of 
cross section increments. 

If KCTYP = 1, skip card 16. 

16 7E10.6 HI 

H2 

BF 
HP1 
HP2 

17* 1415 MATCHG(I) 


Buckling dimension in first direction; zero only 
if no buckling loss is considered. 

Buckling dimension in second direction; zero 
except for one-dimensional slabs. 

Buckling factor (»V): ,/V) for plane bound- 
aries; 2(2. 405/i/3 ) for cylindrical boundaries. 

Perturbed buckling dimension (1st direction); 
zero if HI is not to be perturbed. 

Perturbed buckling dimension (2nd direction); 
zero if H2 is not to be perturbed. 

= 0 Particular material in map is not to be per- 
turbed. 

= material number In compacted sequence be- 
ginning with 1 if it is to be perturbed. 

There will be NMID entries of which NMAT will 
be nonzero running from 1 to NMAT. 


I 


Card 


Format 


Variable 


Description 


Cards 18 to 20 form the material identification map. For convenience, the TDSN map 
may be used but it is not required. 

jjc 

18 7110 NMRA(I) Number of mesh intervals per zone in first direc- 

tion. NZONE1 values. 

jjc 

19 7110 NMZA(I) Number of mesh intervals per zone in second 

direction if NZONE2 > 0. NZONE2 values. 

20 1415 IDM(IJ) Material identification number to include in each 

zone. N ZONE 2 sets of cards (1 if NZONE2 = 0) 
with NZONEl values per card. IDM = 0 if no 
cross sections are read in for the zone (i. e. , no 
perturbation). However, through MATCHG 
(card 17) the map from TDSN can be used here 
without having to zero any IDM entry). 

21* Binary XN Use if KFLUX = 1; real fluxes (NIJ values for 

group 1, then NIJ values for group 2, etc. ). 

5(15, E10. 6) NM, XN Use if KFLUX = 2; NM is the number of mesh 

intervals with the flux level XN. The NM entry 
completing the NIJ total must be <0. 

7E10. 6 XN Use if KFLUX = 3; fluxes for all intervals for 

group 1, then for group 2, etc. 

22* XNA Adjoint flux. Same format as on card 21. 

If KAPROX > 0, read in currents on cards 23 and 24; if NZ > 1, read in cards 25 
and 26 as well. The format must be the same as for card 21. The adjoint currents may 
require reversal, just as adjoint fluxes. 

jjc 

23 XJ First direction real current. 

24* XJA First direction adjoint current. 

Read cards 25 and 26 only if KAPROX > 0 and NZ > 1. 

jjc 

25 YJ Second direction real current. 

26* YJA Second direction adjoint current. 

Read in card 27 only if KPAP = -1. 

27* 7F10.8 PNUF(IG) Use if KPAP = -1. Read in perturbed t'-fission 

cross sections (to be combined with the unper- 
turbed fluxes to obtain the real production F(IJ)). 

28 
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Card 

Format 

Variable 

Description 

* 

27 

(Cont. ) 


See equation (18). NMAT card sets, each with 
NG entries. 

* 

28 

Binary 

F (IJ) 

Use if KPAP ^ 0. Read in the real production, 
either perturbed or unperturbed (see Input 


Notes, p. 29). Binary format from TDSN. NIJ 
values. (If KBEFF = 1 or if KPAP = -1 this 
must be the unperturbed production. ) 

29 Binary FA(IJ) Unperturbed adjoint production. Binary format 

from TDSN. NIJ values. 

Read in card 30 only if KPAP = 1. 

NIJ 

30 E10.6 PAP F(IJ)*FA(IJ)*V(IJ) as given from a previous 

IJ 

PERTRAN problem. 

Read in card 31 only if KD = -1. 

31 7110 KRBC = 0 Not a cell condition; that is, no return cur- 

rent across the right boundary (I = NR). 

= 1 Perfect reflection exists across the right 
boundary (I = NR). 

KLBC Same options for left boundary (1 = 1). 

KTBC Same options for top boundary (J = NZ). 

KBBC Same options for bottom boundary (J = 1). 

Read in cards 32 and 33 only if KBEFF = 1. 

32 7E10.6 BETA(I) Delayed neutron fractions; NDSP entries. 

33* 7E10. 6 DELSP(J) Delayed spectra; NDSP sets of cards, each set 

with (IGLO-IGHI+1) entries (high to low energy). 

Input Notes 

Cross sections are required in the TDSN (ref. 3) format. The TDSN cross sections 
for each group are absorption, j/- fission, transport if P-0 (total if P-0 of P-1), up- 
scattering into the group, within-group scattering, and down- scattering into the group 
in a 7E10. 6 format. The P-1 cross sections (P-1 of P-1) do not occupy the first two 
fields; the rest of the fields contain total, up- scattering, within-group scattering, and 

29 
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down-scattering cross sections. The P-1 scattering cross sections already contain the 
21 + 1 multiplier. The removal cross sections from all energy groups are listed con- 
secutively, 7 to a card. The removal cross sections and the TDSN group sets are 
ordered from high energy to low. 

Within PERTRAN the format is changed to a slight modification of the TDSN format. 
Perturbed cross sections may be used and the increments then calculated within 
PERTRAN. The P-1 cross sections are not considered separate materials and must be 
provided immediately after the appropriate P-0 cross sections for each of the NMAT 
materials (see card 15). 

The real and adjoint fluxes and currents should be obtained from TDSN transport 
calculations - the real and adjoint solutions having been converged to the same multipli- 
cation factor k = k + . All fluxes and currents from TDSN are punched in binary form 
continuously for all mesh intervals for each group. 

The identification map for a perturbation problem may be identical to that used in 
TDSN. PERTRAN thus requires information as to which material regions in this map 
are to be perturbed. 

The normalization factor determined by equation (18) requires the perturbed 
z'-fission cross sections and the unperturbed real fluxes. The KPAP = -1 options forms 


the production quantity Zfa) P <p g from this information. However, the TDSN pro- 

g 

gram provides in convenient form the quantities ^ <p ^ and £K) P v l from 


the unperturbed and perturbed calculations. Under certain conditions these quantities, 
which can be more conveniently handled, may be satisfactory. 

For example, if z'S, is not perturbed or if all vH, are perturbed by the same 

g g 

factor, then the unperturbed production may be used (KPAP = 0). In the latter case this 
factor must then be removed from the printed perturbation results. On the other hand, 

if all the yS f are not perturbed by the same factor but the perturbed fluxes (p^. are 

g g 

not significantly different from the unperturbed fluxes, then the perturbed production 

may be used (KPAP = 0). 

The lifetime is calculated as an absorption perturbation (KP=KS=KD=0; KAR=2) in 
which (1/v) cross sections are supplied rather than absorption cross sections (KCTYP 
may be 0 or 1). Furthermore, the unperturbed real production rate should be used and 
any approximation (KAPROX) may be used. However, if KAPROX = 0, then KCTYP 
must be 1 in order to provide transport cross sections for the Fick's law current 
approximation. Any lifetime "perturbation” should extend over the whole region that 
produced the particular flux spectrum and averaged 1/v cross sections. The fission 
spectrum and the P-1 and removal cross sections are not used so that blank cards may 
be read in for them. 
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The calculation of beta effective requires that the unperturbed production be used 
both separately and combined with the adjoint production; hence, the KPAP=0 option 
must be used. Furthermore, only the adjoint fluxes are needed so KAPROX=0 should 
also be used. In this case KCTYP may be either 0 or 1; no current approximation is 
made. Each of the variables, NMAT, KP, KS, KAR, KD, NZONE1, NZONE2, and 
NMID should be equal to 0 when KBEFF=1. 

Within TDSN the normalization of the fission source in an adjoint calculation results 
in an adjoint flux containing a factor of k. If the prompt fission spectrum is treated as 
the delayed spectrum and a beta effective calculation is performed, then the unadjusted 
adjoint flux provides a /3 g ^ equal to k. Within PERTRAN this factor of k is removed 
(in a calculation) from the adjoint flux so that an absolute is obtained. If a 

delayed neutron fraction of 1. 0 is read in, the ratio /3 g ff//3 is obtained. 

Because it is difficult to determine explicitly from the perturbed and unper- 

turbed forms of equation (25), the calculation of the buckling loss requires an unper- 
turbed Sj. r as well as the increment AS^ r ; thus, the KCTYP=0 option must be used for 
a buckling loss calculation. The perturbation of a buckling dimension H may be treated 
separately or in combination with a perturbation of 2^. 


Output Notes 

The input parameters are listed and the computer storage required for the problem 
is listed under LAST. The incremental cross sections are labeled and listed in the in- 
ternal modified TDSN form. 

If LISTFX = 1 , the regular and adjoint flux, current, and production are given in the 
output. The normalization factor in equation (18) is printed out. 

The contribution to A(l/k) is given by group for production and inscattering sources 
and for absorption, outscatter (removal), and leakage losses. Subtotals provide the 
total contribution by group and type of process. Each GROUP entry is the sum of all the 
listed contributions, with one exception. In the P-1 approximation (KAPROX=2) the 
LEAKAGE contribution has already been included in other categories (see eq. (24)). 

The eigenvalue increment resulting from a change in the buckling loss is listed as 
TRANSVERSE LEAKAGE in the output. The nontransverse leakage out of a reactor sys- 
tem is given by the current- weighted contribution to A( 1/k). For the P-1 approxima- 

5jC 

tion it is listed as LEAKAGE in the output; for the P-0 approximation it is listed as 
TRANSPORT; and for the diffusion approximation it is the sum of the two DIFFUSION 
COEFFICIENT entries. 

The perturbed and the unperturbed multiplications factors (k P and k) are also pro- 
vided. If the lifetime has been calculated it appears as the absorption contribution to 
A (1/k). If a /3 e £j calculation has been performed the delayed spectra are provided as 
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output, and for each delayed spectrum j considered, the corresponding spectral sums 
/3j and /3- ( eff ) are given. 

Programming Notes 

Basically PERTRAN provides one large storage array (X) with a length of 25 000. 
Within this array the number of groups, types of cross sections, materials, and mesh 
intervals are variable. (The amount of this 25 000 storage actually used is listed in 
the output as LAST. ) The dimension of this one large array could be easily changed to 
accommodate modifications to the program. 

Three other easily changed storage constraints occur in PERTRAN: (1) the PNUF 
array (in subroutine INPUT), which must contain (NG*NMAT) locations, is set at 50; 

(2) the STOR array (in subroutine CREAD) must contain the larger of (NTYPS+2* NG) or 
(3* NG) locations and is presently set at 70; (3) the MATCHG array (in subroutine 
IDACAV) must contain NMID locations and is now set at 25. For further convenience in 
modification, an extra common block CALL3 (which is not used in the present version of 
PERTRAN) has been included. 
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APPENDIX D 


PROGRAM LISTING 


This appendix contains the listing (in FORTRAN IV) of the program, overlay infor- 
mation, and a memory map. 


iIBHC PERSN 


THfc C 0 MM J N STATEMENTS 


COMMLJiM 

X 





COMMON 

/CALLI/ 



NG, 

NTYPS , 

NTYNG, 

NF AST, 

NUP , 

NR , 

NZ , 

NI J, 

> 

KF LUX , 

KbEFF , 

KCTYP, 

KIN, 


i 

KP , 

K S , 

KAR , 

KD, 


> KAPRJX,NMAT 

, N M I 0 ,NOSP , IGHI , IGLO 

, X l R , X IZ 



COMMON 

/C ALL2 / 

PAP 




COMMON 

/CALL3/ 

KXTRA1 , 

KXTRA2 , 

KXTRA3, 

KXTRA4 

KXTRA3, 

EXTAA1, 

EXTRA2 , 

EXTRA3 , 

EXTRA4, 

EXT RA 3 

COMMON 

/LINPT/ 

LMA , 

LV, 

LOE LR, 

LOELZ , 



LC , 

LC H i , 

LF, 

LN, 

! LSI A, LXJ, 

L X J A , L Y J, 

LYJA, 

LP, 

LS , 

LA, 

> Lt< , 

LO 1 , 

LO 2 , 

L T LI , 

LTL2, 

LS IX, 

i LR 1 f 

LZ 1, 

LC 1 , 

LF A, 

LNM, 

LAM, 

i L I 0 1 , 

L 102, 

Li OM, 

XKEFF , 



> LAST1, 

LA ST2 , 

LAST 





THE 0 I MEN SIJN STATtMENTS 
DIMENSION X ( 2 5000) 


1000 CALL INPUT 

IF ( L M A ) 1000,1000,1025 


1025 CALL PcKTUR ( X(LMA), X(LV), 
IX KE FF , X ( L N ) , X(LNA) , 

2 X( LD1) , X( LD2) , X( LXJ) , 

3 X(L TL1 ) ,X( LTL2) , X(LS1X)) 


X(LDELR), X(LOELZ), X(LC), XILCHI), 
X ( LF ) , X ( L P ) , X ( LS ) , X ( L A ) , X(LR>, 

X( LXJ A ) , X( LYJ) , X< LYJA) , 


IF( K.APROX .GT.O) GO TO 1050 

IF ( KD .to. O.AND. KAR.NE . 2 ) GO TO 1050 

CALL JIPUSE ( XiLMA) ,X(LV) ,X(LDELR) .X(LDELZ) , X I LC ) , X ( LN I , X ( LNA), 
1 X { LO 1 ),X ( L02) , XILTL1I, XILTL2), X(LA) ) 
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1050 CALL OUTPUT { X(LP), X(LS), X(LA), X(LR), X(LD1), X(LD2) f 
IX ( L YJ A ) » XKEFF t X(LTLi), X(LTL2)) 

GO TO 100 0 
C 
C 

c 

END 


$ I Bl- TC PERI 

SUBROUTINE INPUT 
C 
C 
C 

c 

c 

C THE COMMON STATEMENTS 


CUMMUN 

X 





COMMON 

/CALL 1 / 



NG, 

NTYPS , 

1 NT Y.NGf 

NF AST, 

NUP , 

NR, 

NZ , 

NI J, 

2 

KF LUX , 

KBt F F , 

KCTYP, 

KIN, 


3 

KP , 

K S , 

KAR, 

KD, 


4 KAPRO X f N MAT 

,NMID , NO SP , I GH I , I GLO 

fXIRfXIZ 



COMMON 

/CALL2/ 

PAP 




CUM MON 

/C aL L3 / 

KXTKAi , 

KXTRA2 , 

KXTRA3, 

KXTRA4 

1 KXTRA5, 

EXTRA1 , 

EXTRA2 , 

EXTRA3, 

EXTRA4, 

EXTRA 5 

COMMON 

/CINPT/ 

LM A , 

LV, 

LDELR, 

LDELZ , 

1 


LC , 

LC H 1 , 

LF, 

LN, 

2 LN A * LXJ, 

LXJA,LYJ, 

LYJA , 

LP, 

LS, 

LA, 

3 LR , 

LD1, 

LD2 , 

l r li , 

LTL2, 

LS IX, 

4 LR 1 » 

LZ1, 

LC 1 , 

LF A , 

LNM, 

LRM, 

5 L I 0 If 

L I D 2 , 

L1DM, 

XKEFF , 



6 LASI1, 

LA ST2 , 

LAST 




COMMON 

/CHANG/ 


NO OWN, 



1 NG.M1, 

NGP1 , 


N 




C 

C THE DIMENSION STATEMENTS 

DIMENSION X ( 25000) 

DIMENSION T I TLE< 14) 

DIMENSION PNUF { 50 ) 

C IF KPAP.LT.O, PNUF IS USED AND MUST CONTAIN NG*NMAT LOCATIONS 

C 

C THE FORMAT STATEMENTS 

100 FORMAT ( JLH1) 

10 1 FORMAT! 7F 10. 8) 

102 FORMAT ( F i 0. 6 ? 31 10) 

101 FORMA T { 1HL t 78H THE PERTURBED NU*FISSIGN CROSS SECTIONS ARE (BY GRO 
1UP AND PERTURBED MATERIAL)) 

10 4 FORMA T ( BE 16. 7) 

105 FORMAT { I 1,1X ,14A5) 

106 FORMAT (2X,14A5) 

110 FORMAT (7110) 

111 FORMA T ( 71 16) 
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no no noon 


112 F0RMAT(3I 16, F 16.8, 31 16) 

12 0 FORMAT ( 1HJ, 1 IX , 4HNMA T , 14X ,2HNG , 11 X ,5HNF AST , 1 IX , 5H NDOWN , 13X , 3HNUP, 
1 14 X * 2HNR »14X»2HNZ) 

121 FORMAT (1HJ,10X,5HKCTYP,13X,3HKI N, 1 2X ,4HKPAP , 14X , 2 HKP , 14X, 2HK.S » 

1 13X,3HKAR , 14X,2HtO ) 

122 FORMAT ( 1H J , 9X , 6H NZONE 1 , 1 OX , 6HNZ0NE2 ,1 OX , 6HLI S T FX , 10X , 6H NMID, 10K 
It 6H KAPROX , 10X » 6H KFLUX,10X,6H KBEFF) 

130 FORMAT ( 1 HL , 3HNR= , I 2 , 2 OH IS LESS THAN THREE.) 

131 FORMAT (1HL,14HNZ EQUALS TWO.) 

135 FORMAT { 1 HJ , 1 IX ,4HLAS T ,1 1 X , 5H LAST1 , 1 1 X t5HL AS T 2 , 1 IX , 5HX KEFF , 12X , 
14HNDSP » 12 X t 4H IGHI t!2X»4HIGL0) 


1000 WRI TE (6, 100) 

1005 READ ( 5t 1 05) I TEMP , ( T 1 TLE ( I ) , 1=1,14) 

WRITE (6,106) (TITLE(i), 1=1,14) 

IF (ITEMP) 1005,1005,1010 

1010 REAO (5,110) NMA T , NO, NFAST , NO OWN ,NUP,NR,NZ 
WRITE (6,120) 

WRITE (6,111) NMA T »NG , NF A ST , NDCJWN , NUP , NR , NZ 
REAO (5,110) KCTYP»KIN»KPAP,KP,KS»KAR»KD 
WRITE (6,121) 

WRITE (6,111) KCTYP , K I N,KPAP, KP , KS ,KAR»KO 

REAO (5,110) NZ0NE1 ,NZ0NE2 , LI STF X,NMI 0 , KAPROX , KFLJX , KBEFF 
WRITE (6,122) 

WRITE (6,111) NZJNEl,NZUNt2 , L I S TFX , NMI D , KAPROX , KFLUX, KBEFF 
REAOi 5, 102) XKEFF ,NOSP, IOHI , IGLO 


1050 IF (NZ) 1055,1055,1060 
1055 N Z = 1 

1060 IF (KD) 1085,1085,1065 
1065 IF ( NR — 3) 1070,1075,1 075 

1070 wR l TE ( 6, 130) NR 
RETURN 

1075 IF ( N Z— 2 ) 1085,1080,1085 

1U80 WR I TE (6,131) 

RETURN 

108 5 N IJ =N R *N Z 

NTYPS=NDJ WN+NUP+4 
NTYNG=NTYPS*NG 
NGP 1 = N G <• 1 
NGM 1=NG-1 

IBSTOR =NI J+ND SP* ( IGLJ- IGHU- 1) + NDSP 


1100 LMA = 1 

L V = LMA +NI J 
LDELR = LV+N I J 
L0ELZ=LDELR+NR-1 
IF (NZ-1) 1105,1105,1110 

1105 LC=LDELZ+1 
GO TO 1115 
1110 LC=LDELZ+NZ-1 
1115 lP1CS=LC 

IF( KAPROX. EJ. 2) LP1CS=LC+NMAT*NTYNG 
L CH I = L P 1C S+NM A T*N TYNG 
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C BUCKLING LOSS STORAGE (BELOW] 

LCH1 =LCHI +3*NG*NMAT 
LF=LC H I+NG 
LN=LF +N IJ 
LNA=LN +NG*Nl J 
LXJ =LNA+NG*NI J 
I AOD= 1 

IF( KAPRUX.GT. 0) I A00=NG*N I J 
L XJ A=L XJ+ IAOl) 

L V J =L X JA+ I AOU 
IF(NZ.LE.l) I ADD = l 
LYJA=L YJMAOD 
LP =L Y J A + I AUU 

1 Ft KBEFF.GT.O) lp=lyja+ibstok 

L S = LP +NG 

L A = L S +-NG 

LR = LA tNo 

L 0 1 =L R +NG 

L L) 2 =L Q H-NG 

LTL i=L 0 2+NG 

L TL 2=L TL L +NG* NMAT 

L S i X=L TL2+NG*NMAT 

LAST2=LS1 X+NG- 1 

LAST! = STORAGE REQUIRED BY OVERLAY Oil 
LAST2 = STORAGE REQUIRED BY OVERLAY 001 

1125 LNM=LP 
LRM=LNM+5 
L R 1 =L RM+5 

GO TO ( 1126,1 127,1126) , KIN 

1126 L Z 1=L R 1+NR + 1 
LAST 1=LZ1+NZ 
GO TO 112 8 

1 127 LZ l=LR 1 + 1 
LAST1=LZ1+1 

1128 LC 1 =L P 
ITEMP=LC1+NTYNG— 1 

IF ( I TtMP-LASTl ) 1135 ,1135,1130 
1130 LAST1 = ITEMP 

1135 L F A =L R 1 

IF (KPAP) 1136,1136,1145 

1136 ITEMP=LFA+NIJ-1 

IF ( I TEMP-LAST1) 1145,1145,1140 
1140 L A S T1 = I Tt MP 
1145 L ID 1 = L P 

L ID 2=L ID 1 +NZQNE 1 
L IDM=L I02+NZUNE2 
ITEMP=LIQM+NZ0NE1 
IF < 1 TEMP-LAST1) 1175,1175,1150 
1150 L A S T 1 = I TE MP 


1175 L AST=L AST 2 

IF ( LA ST 1-LAST ) 1185,1185,1180 
1180 LA$T=LAST1 

1185 WRITE (6, 135) 

1186 WRITE (6,112) LAST » LA ST 1 , LAST2 , XKE FF , NDSP , I GH I , I GLQ 

IF (LAST-25000) 1200,1200,1190 
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119G LMA =0 

GO TO 140 0 


1200 CALL VRZ ( X(LNM), X(LRM), XILV), X(LRl), XILZl), X(LDELR), 

1 X(LDELZ) ) 

C 

1225 I Ft KbEPF.EO. 1) GO TU 1275 

CALL CREAD l X(LCHI), XlLC), XILC1) ) 

C 

1250 CALL IOACAV ( (MR, NZ, NZ0NE1 , NZON E2 , NMI D, X l L 1 01 ) , X ( L I 02 > , X ( L I DM ) 
1,X( LMA) > 

C 

1275 CALL iMREAD ( LISTPX, XILNMJ, X(LRM), XILNI, X(LNA), X(LXJ), 

1 X(LXJA), X(LYJ), X ( LYJA) , X(LV) ) 

C 

IP (K.PAP.GE.O) GO TO 1300 
I GO =0 

WR I TE ( 6, 1 03) 

00 230 C I N=i , NMA T 
I SP = I GO+N G 

IGO =N G*< 1N-D + 1 

R tAOl 5,101) (PNUP ( IG) , 1G= IGO, 1 SP) 

2300 WR I TE ( 6, 1 04) ( PN JF ( IG ) , I G = 1 GO , I S P) 

LFENO=LP+NI J- 1 
DO 22 C C I J =LF , LFt ND 
2200 X( I J ) =C.O 

1 JG=LN-1 

DO 200 C 1 G=1,NG 
DO 200 C IJ=1,NIJ 
I J G— I J G+ 1 
K =X ( I J ) 

I F ( K.cU.O ) GO TO 2000 
LA=LP+ IJ- 1 
Ld=NG* (K- 1 ) + I G 

X(LA)=X(LA)+ XI IJGI*PNUP( Lb ) 

2000 CON TIM OE 

C STORING IN FISSION POWER ARRAY F OF SUBROUTINE PAPCAL 

1300 CALL PAPCAL ( KPAP f LISTFX, X ( LNM) ? X ( LHM) » X(LV) f X(LF), 

1 X(LFA) ) 

C 

C 

C 

1400 RETURN 
C 

END 
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$ I8FTC PER2 


SUBROUTINE VRZ 4 NM, RM , V, R, Z, DELR, OELZ ) 


THE COMMON STATEMENTS 




COMMON 

/CALL l / 


NG, 

NTYPS 

I NTYNGf 

NFAST, 

NUP, NR, 

NZ, 

NIJ, 

2 

KFLUX f 

KBEFF , KCTYP, 

KIN, 


3 

KP, 

KS, KAR, 

KD, 


4 K APR OX* N MA F f 

NMI D t NDSP 

, IoHI , 1 GLO , X IR,XIZ 



THE DIMENSION 

STATEMENTS 



DIMENSION 

NM { 1 ) $ 

RM(i) 



0 IMEN SION 

VC I) ? 

Rll) , 

Z< 1 ) 


DIMENSION 

OE LR ( L) , 

OELZ(i) 




THE FORMAT STATEMENTS 
ilO FORMAT (7110) 

115 FORMAT < BE 16. 7) 

120 FORMAT ( 1 HL » 7H0E LR ( I ) ) 

121 FORMAT ( 1HL, 7HDELZ( J) ) 

122 FORMAT ( 1HL . 32HINCORREC T NUMBER OF MESH POINTS= , 15 , 52H INCLUDED TO 
1 OBTAIN EITHER R OR Z. SHOULD HAVE BEEN* 1 5 ) 

130 FORMAT ( 1 HL , 5HKGE 0= , I 3 ) 

131 FORMAT ( 1 HL * 6H V ( l » J ) ) 

132 FORMAT ( 1HL,4HR (I ) ) 

133 FORMAT ( 1HL * 4HZ ( J ) ) 

THE FUNCTION STATEMENTS 

I NO EX ( LENGTH, I NOE XO, I NOE XL) =LENGTH * { I NDEXG-1 ) + I NDEX L 


1000 NRP 1 = N R + 1 
NZP 1=NZ+1 

GO TO (1005,1005,1050), KIN 

1005 CALL OCREAD ( V(l), V(NIJ)) 

IF(KBEFF.EU.l) GO TO 1100 
GO TO ( 1006,1100,1050), KIN 

1006 CALL BCREAO < R ( 1 > , R ( NR Pi ) ) 

IF (NZ-2) 1100,1100,1010 

1010 CALL 3CRE AD (Z(l) ,Z(NZP1) ) 

GO TO 1100 

1050 R EAD (5,110 KbEO 
WRITE (6,130) KGEO 

CALL RZREAD ( NM, RM, ITEMP, R ) 
IF (ITEMP-NR) 1055,1060,1055 
1055 WRITE (6, 1221 ITEMP, NR 
GO T 0 120 0 

1060 IF (NZ-1) 1075,1075,1065 

1065 CALL RZREAD ( NM, RM, ITEMP, Z ) 
IF ( I TEMP— NZ ) 10/0,1080,1070 

1070 WRITE (6,122) ITEMP, NZ 

GO TO 120 0 
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1075 TEMP-1.0 
J=1 

GO TO 108 6 
1080 J = 1 

1085 TEM P=Z (J+1)-Z(J) 

1086 £J = INDfcX( NR, J ,0) 

GO TO (1087,1088,1089), KGhU 

1087 AVE=1.C 

GO TO 1090 

1088 A VE = 3 . 141 59265 
GO TO 109 0 

1089 A Vt = 3. 141 59265*4. 0/3, 0 

1090 DO 1091 1 — 1 , NR 
IJ=IJ + 1 

1091 VC 1 J ) = fEMP*(R ( I + i)**KGEG-R( I ) **KG£U) 
IF (NZ-J) 1100,1100,1092 

1092 J=J+1 

GO TO 108 5 
C 

c 

1100 WRI TE (6, 131) 

DO 1105 J = 1 , N Z 
I = INDEX(iMR, J f 1) 

I TEMP = 1NDE X ( NR , J , NR ) 

1105 WRITE (6,115) (V(IJ), IJ=I,ITEMP) 

IF( KBEFF.tO. 1 ) GO TO 1200 
GO TO (1110,1150,1110), KIN 
1110 WRI TE (6, 132) 

WRITE (6,115) (R(I), 1=1, NRPi ) 

XIR=R( 2)/2. 

IF (NZ-2) 1150,1115,1120 

1115 GU TO ( 1150, i 150,1120) , KIN 
1120 WRITE (6, 133) 

WRITE (6,115) (Z(J), J =1 » NZP1 ) 

X I Z =Z ( 2 )/ 2 . 

C 

1150 I TEMP = NR— 1 

GO TO ( 11 55, 1165,1155), KIN 
1155 DO 11 6 C I =1, I TEMP 
1160 DEL R( l )=0 * 5* ( R ( 1 + 2 ) -R ( 1 ) ) 

GO TO 1170 

1165 CALL dCREAQ ( DE LR ( L ) , DE LR ( 1 TE MP ) ) 
1170 WRI TE (6, 120) 

WRITE (6,115) (DELR(l), I = 1,ITEMP) 

IF (NZ-2) L200, 1200, 1175 

1175 I TEMP = NZ- 1 

GO TO (1180,1190,1180), KIN 
1180 00 1185 I =1, I TEMP 
1185 DELZI I )=0.5*< Z(I+2)-Z(l)) 

GO TO 119 5 

1190 CALL 6CREA0 ( DE LZ { 1 ) , DE LZ ( I TE MP ) ) 
1195 WRI TE (6, 121 ) 

WRITE (b,115) (DELZ(I), I=1,ITEMP) 

C 

c 

c 

1200 RETURN 
C 

END 
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*IBFTC PER21 


SUBROUTINE RZREAD ( N M, KM, 
C 
C 
C 
C 

c 

C THE DIMENSION STATEMENTS 

DIMENSION NM(1), 

DIMENSION Rd) 

C 

C THE FORMAT STATEMENTS 

116 FORMAT ( 5( It), E10.6) ) 

C 

C 

C 

100C R ( i) = 0.0 
K OUN T = 0 
K STOP = 1 

1005 READ (5,116) (NM(IJ f RM(I), 
DO 1025 1=1, 5 
IF (NM( I) ) 1010,1025,1015 

1010 NM{ il=-NM(I> 

K STOP = 2 

1015 K 1 =K0 UN T+ 1 

KGUN T = KOUNT+NM ( I ) 

0 EEL =R M ( I )-R(Kl) 

T EM P=NN ( I ) 

DEL L = 0 ELL / TE MP 
DU 1020 K =K 1 , ROUN T 

1020 R(K + 1)=R< Kj+DELL 

GO TO (1025,1050), KSTGP 

1025 CON T IN LE 

GO TO 100 5 
C 
C 
C 

1050 RETURN 
C 

END 


KOUNT, R ) 


RM( 1) 


1-1*5) 
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$ IBFTC PER3 


SUBROUTINE CREAD ( CHI, C, Cl » 
C 
C 

C READ IN CROSS SECTIONS 

C 

C 

C THE COMMON STATEMENTS 


COMMON 

/CALL1/ 



NG f 

NTYPS, 

1 NTYNG, 

NPAST, 

NUP * 

NK, 

Hit 

N I J, 

2 

KFLUX, 

KBtFF , 

KCTYP, 

KIN, 


3 

KP , 

KS » 

KAK, 

KD, 



4 KAPROX,NMA T, NMID ,NOSP , 1GHI , I GLU ,X I R , X IZ 
COMMON /CHANG/ NDOWN, 

1 NGM l , NGP1, N 

THE DIMENSION STATEMENTS 

DIMENSION C(I), CHI (II, Cl (1) 

DIMENSION STOR (70) 

STJR MLST CJNTAIN THE LARoER OF ( NT YPS + 2 *NG) OR <3*NG) LOCATIONS 

THE FORMAT STATEMENTS 

101 FORMAT ( 1H ) 

102 FORMAT! 1HK.36H DELTA BUCKLING LOSS CROSS SECTIONS ) 

ID 3 FORMAT! IHJ, 5X ,1 BH FIRST DIRECTION /(7E18.7)) 

104 FORMAT! 1HJ*5X*18H SECOND DIRECTION /(7E18.7I) 

105 FORMAT! 1HK,40X,10H MATERIAL ,121 

106 FORMAT! 1HR, 31H BUCKLING LOSS INFORMATION ****/18H FIRST DIMENSION 

1 , F ID , 6 , 1 8H StCOND DIMENSION ,F10.6/2dH FIRST DIMENSION PERTJRBED 

2 , F 10 . 6, 2 8H SECOND D 1 Mt NS I ON PERTURBED .F10.6/18H BUCKLING FACTOR 

3 , F 12 . 7 I 

112 FORMAT ( TtlO.o) 

113 FORMAT ( IE Id. 7) 

120 FORMAT ( 1 HL , 7HCH l ( I G I I 

121 FORMAT ( 1 HL , 2 6HRP 1 FORMAT CROSS SECTIONS.) 

122 FORMAT ( 1HL , 1 2 5HC IN , l G , I T Y ) FOR ITY OF S1G A, NO SIG F, ( 0* OR TR 3 
1R TUT), SIG OUTSCATTER, NOP VALUES OF SIG UP, AND NDOWN VALUES OF 
2SIG DUwN. ) 

123 FORMAT ( 1 HL , 8 OH TD SN FORMAT CROSS SECTIONS. TRANSPORT APPROX IMAT ION 
1. P-0 TRANSPORT CORRECTED CS.) 

124 FORMAT (1HJ,80HTJSN FORMAT CROSS SECTIONS WITH D* IN PLACE OF SIG 

1 TR . J IFF U SI UN APPRUXI MAT1 UN. /83H U*= ( DELT A S 1GTR/ ( 3.*! S1GTR*<=2) ) ) 

2 OR D *= ! - DEL TA DIFCOEF/! l.+DELTA Di FCOEF/DI FCOEF) ) i 

125 FORMAT! 1HL ,91HDEL TAS COMPUTED INTERNALLY FROM DIFFERENCE BETWEEN P 
1ERTUR3ED AND UNPERTURBED CRUSS SECTIONS) 

126 FORMAT ( IHJ , 73HTDSN FORMAT CROSS SECTIONS. TRANSPORT APPROXIMATION 
l. WITH P-0 AND P-1 C S. ) 

127 FORMA T ( 1HL , 1 5HDELTA S READ IN ) 

128 FORMAT (1HK,59H TD SN FORMAT. P-1 CROSS SECTIONS WITH FACTOR OF 3 IN 
1CLUDED ) 

129 FORMAT! 1HK,33H DE L TA OUTSCATTER CROSS SECTIONS ) 

130 FORMAT ( 1HL , 1 1 6H C(N,IG,ITY) FOR ITY OF 0.0 , P-0 G-TO-G, TRANSPORT 

1, P-1 G-TU-G, NUP VALUES OF SIG UP, AND NDOWN VALUES OF SIG DOWN.) 

THE FUNCTION STATEMENTS 
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INDEX ( LENGTH* INDEXO,! NOE XL) =LENGTH*( l NDEXO-1 ) + I NDE X L 
C 

c 

c 

1000 READ ( 5,112) (CHIUG), IG=1,NG) 

WRITE (6,120) 

WRITE (6,113) (CHKIG), IG=1,NG> 

NPOCS=NTYNG*NMA T 
C JQG=0 P-0 CROSS SECTIONS 

C J0G = 1 P-1 CROSS SECTIONS 

DO 290 H- 1 , NMA T 
NO WP 1 = C 
JQG=0 

I F ( KC T YP • NE • 0 ) GO TO 1025 

279 JDOUBL-1 

C JDOUBL = 1 READ IN PERTURBED CROSS SECTIONS 

C J DOUBL =2 READ IN UNPERTURBED CROSS SECTIONS 

281 DO 280 JIG=1,NG 

JK 1 = INDEX ( N TYP S ? J IG f 1 ) 

J K2= IN DEX ( N TY PS, JIG tNTYPS) 

IF( JDOCBL .EQ. 2) GO TO 282 
READ( 5,112) (Cl(JK) ,JK=JK1, JK2I 
GO TO 280 

282 R EAD( 5,112) ( STOR ( JI K ) , JI K=1 ,NTYPS) 

HULD^STOR ( 3) 

L END = 7 C— N G+J I G 
STUR( LEND ) =ST0R(3) 

KUN T= 0 

DO 270 RK = JK 1 , JK2 
KON T=K JNT + 1 

Cl( KK ) =C 1 ( KK ) - S TOR ( KONT) 

I F ( AB S ( C 1 ( KK ) ) .LT.. 0000005) C1(KK)=0.0 
27 0 CONTINUE 

I F < KAPRGX*GT *0) GO TO 280 

Cl( JR 1+2)= C 1( JKL+2)/(3.0*H0LD*HGLD) 

C SEE EQ. 32 IN WRITEUP 

280 CONTINUE 

IF( JDOCBL .EU. 2) GO TO 210 
JDOUBL =2 
GU TO 281 
C 
C 

1025 GU TO (1031,1100), KCTYP 
210 IF( JOG.EsJ.il GO TO 1032 
J QG = 1 

WRI TE ( 6, 125) 

I F ( KAPRUX.EQ.O) WRI TE (6,124) 

IF( KAPROX.EO. 1) WRI TE (6,123 ) 

I F ( KAPR0X.EQ.2) WRI TE (6,12b) 

GU TO 103 2 

1031 WRITE (6, 127) 

1032 WRI TE ( 6,105) N 
DO 1055 l G=1 , NG 

K 1= IN DEX( N TYP S , IG , 1 ) 

K 2= INDEX! NTYPS, IG,NTYPS) 

IF(KCTYP.EQ.O) GO TO 300 
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READ (5*112) (C 1 ( K) * K=Ki,K2) 

300 WRITE(6»113) (C l ( K) , K=K1 , K2 ) 

TRANSFER TDSN CROSS SECTIONS FROM INPUT ARRAY Cl TO PERMANENT 
STORAGE (ALTERED FORMAT) IN C 
L = INDtX(N TYNG * N, K1 ) 

LL=L+ 1 

IF ( N0WP1.EQ. 1 ) L=L+NP0CS 
DO 1055 K =K 1 » K2 
LLL =NP CCS+LL 

IF1N0WP1.EU. l.AND.K.EQ. (Kl+1) > GO TO 1055 
I F ( K.GT.( Kl + 2) ) GO TO 1045 
1035 C(LKIIK) 

GO TO 105 5 

1045 IF(N0WPI. EQ. l.ANJ.K.EU. (Kl + NUP+3) ) C ( LLL+2 ) = Cl ( K. ) 

I F ( NO WP 1. EQ • 1) GO TO 1046 

I F ( KAPROX.EQ. 2. AND. K. EQ . (Kl+NUP+3) ) C(LLL)=C1(K) 

1046 IF IK- (K1+NUP+3I) 1050,1055,1035 

1050 C ( L + 1 ) =C 1 (K) 

1055 L =L +1 

IFIKAPROX .Lt. 1) GO TO 1054 
I F ( NUWP1.EQ. I ) GO TO 1054 
WR I TE( 6,128) 

NO WP 1=1 
GO TO 279 
1054 WR I TE ( 6, 101) 

READ OCTSCATTER 
WRITEl 6, 129) 

READ (5,112) ( C 1 ( K) , K=1 , NG ) 

I F ( KC TYP.NE.O) GO TO 250 
READl 5,112) ( STOK ( K ) , K=1 , NG ) 

DO 240 JX = 1, NG 
240 Cll JX) =Ci( JX)-STOR( JX) 

250 WRl TE ( 6, 1 13) ( C 1 ( K) , K-=l , NG ) 

DO 10 o C I G = i * NG 
K 1= INDEX( NTYP S , IG,4) 

L = INDE X(N TYNG ,N ,K1 ) 

1060 C(L )=C 1( IG ) 

H 1 = 0 . 0 

IF(KCTYP.EQ.l) GO TO 290 

READ TRANSVERSE LEA KA l,E (8UCKLING) LOSS INFORMATION 
CALCULATE BUCKLING LOSS CROSS SECTIONS AND STORE IN C 
GAMMA = 0.71045608 
I TAG=0 

lFtN-GT-1) GU TO 1200 

READl 5,112) HI ,H2,BF ,HP1 ,HP2 

WR I TE ( 6,1 06) HI ,H2 ,HP l ,HP2 ,BF 

I F ( Hl.EQ. 0.0) GU TO 1150 

H = H 1 

HH = H2 

I ST =7 0— NG 

NS=NG*NMAT*NTYPS 

IF ( KAPRUX .EQ . 2) NS=2*NS 

I A=N S + 1 

I8=3*NG*NMAT+IA-1 
DO 120 1 I =IA, IB 


no non 


I I l I 


i 


min i i ii ii 


II III Mil II III I ■ II II III I II I 1 1 III III II ■ III III I ■IIIHII III HIM 


mu mi 


1201 Cl 1 )=0.0 

1200 DO 1202 I LK=1 ,NG 
I TAG= I TAG + 1 
INDW=I ST+ILK 
IK=NG*(N- 1KNS+ILK 

HANGUN = (H*ST0R.( I NCIW ) + 1 .42 0912 16 ) **2 
Cl IK ) = 1 STUR ( I NCJW)*BF*BF/HANGGN) -Cl IK) 

IF! H2 .EU . 0. 0) GO TO 1203 
HANo0N=l HH*S TORI I NOW) + 1. 42091 216) 

IKK=NG*NMAT+I K 

Cl IKK )=( S TOR 1 IN0W)*BF *BF /HANG ON ) -C (IKK) 

1203 JJ=INDfcX( NTYP S» ILK, 3) 

JJJ=lNDtXlNTYNG,N,Jj) 

IF(KAPR0X.EQ.2) J J J=J J J+NPOC S 
IF(KAPROX.Gt.i) FAC TOR = 1 . 0 

IF! KAPROX.EO. 0) F AC T0R=3. * l STOR 1 INQw) **2) 

IKkK=2*NG*NMAT+IK 

STORE 1/1 3* 1 S1GTR**2) ) IN C — USED FOR CURRENT (FICK.S LAW) 
APPROXIMATION FUR KAPR0X=0. J = -DELPHI / (3*S IGTR) 

C( IKKK ) = 1 .0/FAC TOR 

STUR( I NOW ) = S T OR ( I NOWI+FAC TOR*C( JJJ) 

1202 CONTINUE 

I Ft HP 1 .NE .0.0) H=HP 1 
IF< HP 2.NE .0.0) HH =HP2 
IF( ITAG.EQ.NG) GO TO 1200 
290 WRI TEl 6,101) 

GO TO 1130 

PROVISION FOR DIFFUSION THEORY CKUSS SECTIONS 
1100 WRITE 1 o» 121) 

DO 1105 N =lf NMA T 
CALL XCHANG (Cl. C ) 

1105 WRITE 16, 101) 

C 

1150 WRITE (6, 122) 

J T = 0 

1151 00 1160 N = 1, NMA T 

WR l TE 1 6,105) N 
00 1155 I G=1 , NG 

K 1= INDEX ( N TYP S , l G , 1 ) 

K 2= INO EX (NTYPS,IG,NTYPS) 

K 1= INO EX 1 NT YNG »N » K1 ) 

K 2= INO E X ( N T Y NG , N , K2 ) 

IF(JT.EQ.O) GO TO 1156 
K 1=K 1 +NP0C S 
K2=K2+NP0CS 

1156 WR I TE 1 6, 1 13) 1C ( K) , K=K1 , K2) 

1155 CONTINUE 

IF! KAPROX .ED. 2. AND. JT.EQ.O) GO TO 1160 
I F ( Hl.EU. 0.0) GO TO 1160 
WR I TE 1 6, 102) 

I ST = I A 

1SP=I SI+NG-1 
I1ST=NG*NMAT+IST 
.1 ISP = 1 IST+NG- 1 


44 


non onooooon ooo o non 


WR I TE ( 6, 1 03 ) ( C ( I G) *IG = IST»ISP) 

IF( H2.NE.0.0) WR I TE ( 6 , 1 04 ) ( C ( I G ) , I G= I I ST , 1 1 S P ) 

1160 WR ITE (6, 101) 

IFIKaPROX.LE. 1.0R.JT.EQ.1) GO TO 1161 
WRI TE ( 6, 130) 

JT=1 

GO TO 1151 


1161 RETURN 
END 


16FT 


140 0 


. PER4 

SUBROUTINE XC HANG I Cl, C ) 

PRUVISION FOR DIFFUSION THEURY CROSS SECTIUNS 

RETURN 

END 


IBFTC PEK5 


DECK 


SUBROUTINE IDACAV 
IMAP ) 


SUBROUTINE IDACAV 


( NR , NZ , NZONRA »NZUNZA,NMID, NMRA.NMZA, I DM, 

OBTAINS THE IDENI FICATI ON NUMBER MAP. 


THE DIMENSION STATEMENTS 

DIMENSION NMRA(l), NMZA(l), IDM(l) 

DIMENSION IDMAP(l), MA TCHG ( 25 ) 

MATCHG MUST CONTAIN NM ID LOCATIONS 

THE FORMAT STATEMENTS 
110 FORMAT (7110) 

115 FORMAT (1415) 

116 FORMAT (4313) 

120 FORMAT (1HL.11HMAP MA ( I , J ) I 

121 FORMAT! 1HL.39HMAP MA(I,J> WITH ONLY PERTURBED ENTRIES ) 

126 FORMAT! 1415) 

127 FORMAT ( 1HK,37HTHE NON-ZERO MATERIALS ARE PERTURBED /24I5) 


o o n on on 


C 

c 


THE FUNCTION STATEMENTS 

INDEX ( LENGTH, I NDE XO , I NDE X Li =L£NGTH* ( INDEXO-i ) +INDEXL 


C 

C 

c 

N IJ=NR*NZ 

READ( 5 , 12 6 ) ( MA TCHG { I ) , I =1 , NMI D) 

WRI TE( 6, 12 7) ( MA TCHG ( 1 ) , I =1 , ?NMIU ) 
1000 READ (5,110) ( NMR A ( I ) « I=1,NZ0NRA) 

IF (NZQNZA) 1005,1005,1010 
1005 N ZONZ A = 1 
NMZ A ( 1 ) = 1 
GO TO 102 5 

1010 READ (5,110) (NMZA(J), J=1,NZ0NZA) 

C 

1025 K = 0 
J =0 

DO 105C J J = 1 , NZONZA 

READ (5,115) ( I DM ( I I ) , II=i,NZONRA) 

K=K+NM ZA( JJ) 

1030 J=J+1 

I =0 

DO 10 A 5 I I =1 , N ZUNRA 
L =NMR A (II ) 

DO 1035 K K = 1 » L 
1 = 1+1 

I J = IND EX ( NR # J , l ) 

1035 I D M A P { IJ ) = IOM ( I I ) 

1045 CONTINUE 

IF (J-K) 1030,1050,1050 

1050 CONTINUE 


1075 WRI TE (6, 120) 

J DO G= 0 

J DOG= 0 REGULAR MAP WITH NMID MATERIALS 
J DOG= 1 MAP WITH NMAT PERTURBED MATERIALS 

1076 DO 1030 J = 1, NZ 
K = INDE X{ NR , J , 1) 

L=INDEX(NR, J ,NR ) 

1080 WRITE (6,116) (IDMAPUJ), IJ=K,L) 

I F ( JDO G* t Q * 1 ) GO TO 1100 
DO 12 0 C I K=1,NI J 
KKK = I DMAP ( IK) 

1200 IDMAPI IK) =MATCHG(KKK) 

J DO G= 1 

WRI TE ( 6,121) 

GO TO 1076 


HOC RETURN 
C 

END 
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SIBFTC PER61 


SUBROUTINE READIT ( NMR, RMI , IJGS, NIJ, SV, KENQ I 


THE DIMENSION STATEMENTS 

DIMENSION NMR ( 1 ) , RMI(l), SV(1) 

THE FORMAT STATEMENTS 

118 FORMAT (I 5,E10.6,I5,E10. 6,I5,E10.6,I5»E10.6, 1 5 , E 10 .6 ) 

119 FURMAT ( 1HL » 5 1H TUQ MANY VALUES F CR SUBSCRIPTED VARIABLE SV READ IN 

1 ) 


100C L =0 

KEND=0 

1005 READ (5,118) ( NMR ( I ) , RMI (I) , 1=1,5) 

1 = 1 

1010 IF (NMR(I)) 1050,1050,1015 
1015 L i=L + 1 

L =L +NMRI I ) 

IF (L-NIJ) 1025,1025,1020 
1020 WRI TE (6, 119) 

K EN D= 1 
GU TO 110 0 
1025 I TEMP 1=L1 +1 JGS 
I TEMP 2=L + IJGS 
DU 1030 I JG = I TEMPI, I TEMP2 
1030 SV( IJG )=RMI ( I ) 

IF (L-NIJ) 1035,1100,1100 
1035 IF (1-5) 1040 ,1005,1005 

1040 1=1+1 

GU TO 1010 

1050 IF (L-NIJ) 1055,1100,1100 
1055 IF (L) 10 80,1060, 106 5 
1060 ITEMP i=IJ GS+1 
GO TO 1070 
1065 ITEMP 1=ITEMP2+1 
1070 ITEMP2 = IJ GS+N I J 

DO 1075 I JG=I TEMPI , l TE MP2 
1075 SV( IJG >=0.0 


1100 RETURN 
C 

END 
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$ IBFTC PER6 


SUBROUTINE NREAD { LISTFX, NM, RM, XN, XNA, X J ,X J A, Y J , YJ A, V ) 
C 
C 

c 

c 

c 

C THE COMMON STATEMENTS 


COMMON 

/CALL l / 



NG, 

NTYPS f 

1 NTYNG, 

NF A ST, 

NUP, 

NR « 

N It 

NIJ, 

2 

KFLUX, 

KBtFF , 

KCTYP, 

KIN, 


3 

KP , 

KS , 

KARt 

KOt 



4 KAPRUX,NMAT,NMIO,NDSP, IGHI , l GLO »X IR , X I Z 
C 

C THE DIMENSION STATEMENTS 

DIMENSION NM(i), Vll), RM(1) , XJ<1), XJA(l) 

DIMENSION XN ( 1 ) » XNA(l), YOU), YJA(I) 

C 

C THE FORMAT STATEMENTS 

115 FORMAT ( BE 16. 7) 

116 FORMAT ( 4H IG = , 13) 

117 FORMAT ( 4H J=,I3) 

120 FORMAT ( 1 HL * 1 8H THE FLUXES XN(I,J)) 

121 FORMAT ( 1 HL » 2 7H THE ADJOINT FLUXES XNAd.JH 

122 FORMA T ! 1H L » 3 6H THE FIRST DIRECTION CURRENTS XJ(I,JI) 

123 FORMAT! 1HL,45HTHE FIRST DIRECTION ADJOINT CURRENTS XJA<I,J)) 

124 FORMAT ( 1HL» 37HTHE SECOND DIRECTION CURRENTS YJ(I,J)) 

125 FORMAT (1HL, 46H THE SECOND DIRECTION ADJOINT CURRENTS YJA!l,J)) 

C 

C THE FUNCTION STATEMENTS 

INDEX! LENGTH, I NDE XO , l NDE XL) =LENG TH* ( I NDEXO-1 ) + I NDEXL 
C 
C 

C READ IN THE REGULAR THEN THE ADJOINT 1) FLUX 2) FIRST DIRECTION 

C CURRENT 3) SECOND DIRECTION CURRENT 

C 

IF! KBEFF.EQ.O) GO TO 190 
CALL REED l XNA) 

KR I TE = 2 
GO TO 1106 

190 CALL REEDIXN) 

CALL REED! XNA) 

IF! KAPRUX.EQ. 0) GO TO 1100 
CALL REED(XJ) 

CALL REEDtXJA) 

NIT=NIJ*NG 
DU 200 J= 1 *N I T 
YJ! J)=XJ( J) 

200 Y JA ! J ) =XJ A ( J ) 

IF1NZ.LE.1) GO TO 1100 
CALL REED l YJ ) 

CALL REEDIYJA) 

C 

C 

C 

1100 IF (LISTFX) 1200,1200,1104 

1104 KR I TE = 1 

1105 WRI TE ( 6, 120) 

GO TO 130 0 
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1106 WRITfcC 6» 1 2 1 ) 

GU TO 130 0 

1107 WRI TE ( 6, 122) 

GO TO 130 0 

1108 WRI TE( 6, 123) 

GO TO 130 0 

1109 WRI TE ( 6t 124) 

GO TO 130 0 

1111 WR I TE { 6,125) 

1300 00 111C IG-1,NG 
WRITE (6,116) IG 
DO 1110 J=1,NZ 
WRITE (6,117) J 
K1=INDEX( NR , J , 1) 

K2=INDEX(NR, J,NR) 

K 1= INL) EX ( NIJ, I G , K 1 ) 

K 2= IN 0 EX ( NIJ,IG,K2) 

GO TO < 1201, 1202, 1203 ,1204 ,1205 ,1206 ) , KRITE 

1201 WRI TE ( 6,115) ( XN ( I J ) , I J = K1 , K2 ) 

GO TO 1110 

1202 WRI Tbl 6,1 15) ( XNA ( I J) * I J = Kl , K2 ) 

GO TO 1110 

1203 WR I TE ( 6 , 1 1 5) (XJ ( I J i , 1 J = Kl , K2 ) 

GO TO 1110 

1204 WK I TE ( 6 , 1 1 5 ) (XJA (IJ) , I J- K1 ,K2) 

GO TU 1110 

1205 WRI TE ( 6,115) (YJ (IJ) , I J - K1 , K2 ) 

GO TO 1110 

1206 WR 1 TE ( 6 , 1 1 6 ) (YJA (I J) » I J = K1 ,K2) 

1110 CONTINUE 

KRi TE =KR I TE + 1 
I E( KR I TE . EQ. 2 ) GO TO 1106 
IF{ KAPROX *EQ . 0) GO TO 1200 
I F( KR I TE.EQ. 3 ) GU TO 1107 
IF(KK ITE.EQ.4) GO TO 1108 
IF4NZ .Lfc. 1) GO TO 1200 
IF( KK ITE.EQ. 5) GO TO 1109 
IF{ KR I TE.EQ.6) GO TO 1111 


1200 RETURN 
C 

END 
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SIBFTC PER666 


SUBROUTINE REED ( XYZ ) 
COMMON /CALL 1 / 


NG, 

NTYPS 

NTYNG, NFAST, 

NUP, NR» 

NZ, 

NIJ, 

! KF LUX * 

KBEFF , KCTYP, 

KIN, 


i KP, 

KS, KAR, 

KD, 


» KAPROX,NMAT * NMI0»NDSP 

, IGHI , 1 GLO «X t R t X I Z 




DIMENSION XYZ(l) 

INDEX! LENGTH, I NDE XO , I NOE XL) =LENGTH* ( I NDEXO-L ) + I NDEX L 
100 FORMAT ( 7E 10.61 
C 
C 

GO TO (200,200, 160,150) , KF LUX 
200 DO 100C I G= 1 , NG 

K 1= IND EX( NIJtlGfl) 

K 2= INDEX! NIJ, IG ,NI Jl 
1000 CALL BCREAD!XYZ(K1)»XYZ(K2) ) 

RETURN 

150 DO 200 C I G = 1 » NG 

K 1= IND EX( NIJ, I G , 1 ) 

K 2= INDEX! NIJ, IG ,NI J i 
2000 R £ AD( 5 , 10 0 ) ( XYZ ( I ) ,1 =Ki , K2 ) 

RETURN 
16 C Kl^O 

DO 300 C I G = 1 , NG 

CALL READ ITINMtRMvKlvNIJ rXYZtKE ND) 

IF(KEND) 3000,3000,1700 
300 C K L-K 1 +N I J 
1700 RETURN 
END 


SIBFTC PER7 

SUBROUTINE PAPCAL ( KPAP, LISTFX, NM» RM, V, F, FA ) 


THE COMMON STATEMENTS 


COMMON 

/CALLi/ 


NG, 

1 NTYNG, 

NF A ST, 

NUP, NR, 

NZ, 

2 

KF LUX , 

KBEFF, KCTYP, 

KIN 

3 

KP , 

KS, KAR, 

KD, 

A K APR □ X* N 

MAT,NMID , ND SP 

, IGHI , IGLO,XIR,XIZ 


COMMON 

/CALL 2/ 

PAP 


THE 0 1 MEN 

SION STATEMEN 

TS 


DIMENS ION 

NM ( 1 ) , 

RMIll 


DIMENSION 

V( 1) 



0 IM EN S ION 

F ( 1 ) , 

FA11I 



NTYPS , 
N IJ, 
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THE FORMAT STATEMENTS 
112 FORMAT (7E10.6) 

115 FORMAT (8E16.7) 

117 FORMAT (AH J=,I3) 

120 FDRMAT ( 1HL * 24HTHE POWER DENSITY F(I,J)) 

121 FORMAT ( 1HL i 33HTHE ADJOINT POwER DtNSITY FA( 1 , J ) ) 

15C FORMAT ( 1HL, 32HADJ0I N T POWER TIMES POWER E DUALS , El 6 .8 ) 

THE FUNCTION STATEMENTS 

INDEX! LENGTH, I NDE XO ,1 NDE XL I =LENGTH* ( I NDEXO-1 ) FI NDEXL 


1000 IF (KPAP I 1015*1010,1150 

1015 DO 1016 I 1=1, NIJ 

1016 FA< II )=F( II ) 

1010 CALL BCREAD ( F ( 1 ) * F ( N I J ) ) 

I F ( K.PAP.GE .0) GO TO 1018 

DO 1017 J J = 1 , N I J 

IFCFAI JJ) .NE.0.0) F ( J J ) =F A ( J J ) 

1017 CONTINUE 

1018 CALL BCREAD ( FA C 1 ) ,FA ( NI J )) 


1100 IF (LISTFX) 1 125,1 125,1105 

1105 WRI TE (6, 120) 

DO 11 10 J =1,NZ 
WRITE lb, 117) J 
K 1=IN0EX( NR, J, 1) 

K 2= INDEX! NR,J,NR) 

1110 WRITE (6,115) (F ( I J ) , IJ=Ki,K2) 
WRI TE (6, 121) 

DO 1115 J = l » N Z 
WRITE (6,117) J 
K 1= INO EXl NR , J , 1) 

K2=INDEX(NR, J,NR) 

1115 WRITE (6,115) (FA(IJ), IJ=K1,K2) 


1125 PAP=0 . C 

DO 1 1 3 C 1 J = 1 , N I J 
1130 PAP=PAP+F ( IJ)*FA{ IJ)*V(IJ) 
GO TO 117 5 


1150 READ ( 5,112) PAP 


1175 WRITE (6, 150) PAP 


1200 RETURN 
C 

END 
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SUBROUTINE PERTUR ( MA , V f D£ LR , 0 fc LZ , C t CH I , X KEFF, XN, X NA, F, 
1 P, S, A, R, Ul f D2 , XJ, X JA * YJ, YJA, TLL, TL2 * SIX) 

C 

C TO COMPUTE PERTURBATION SOURCES AND LOSSES * 

C 

c 

c 

c 

c 


THE COMMON 

STATEMENTS 




COMMON 

/CALL 1 / 


NG, 

NTYPS , 

1 NTYNG, 

NFAST, 

NUP, NR, 

NZ, 

NIJ, 

2 

KF LUX , 

KBEFF, KCTYP, 

KI N, 


3 

KP , 

KS, KAR, 

KD, 


4 KAPKOX*M MAT ,NMIU,NOSP 

, IGHl , I GLO ,XIR»XIZ 



COMMON 

/CALL2/ 

PAP 



COMMON 

/CALL3 / 

KXTRA1, KXTRA2, 

KXTRA3, 

KXTRA4, 

1 KXTRAb, 

EXTRA1 , 

EXTRA2, EXTRA3, 

EXT RA4, 

EXTRA5 

THE DIMENSION STATEMENTS 



D IMENSIUN 

MAUI, 

VII) , 



1 

OhLK(l) , 

DELZ (1 ) 



i) IMENSION 

C (1) , 

CHI ( 1) 



DIMENS ION 

XN(l) , 

XNA(l) , 

FI I ) 


DIMENSION 

P ( 1) , 

s m , 

aid, 


1 

R< Li , 

Dim, 

02(1), 

TL1C1), TL2I1) 

DIMENSION 

XJ(1), XJ A ( 1 ) , Y J ( 1 ) , YJA ( 1 ) , 

S1XII ) 



C THE FUNCTION STATEMENTS 

INDEX { LENGTH f I NDE XO , I N DE X L I ^LENGTH* { I NDEXO-1 ) + I NDE X L 
C 
C 

101 FORMA T ( 7E 10.6) 

102 FORMA T { 1H 1 , 5H THE ,12, 40H DELAYED SPECTRA ARE (FROM ENERGY GROUP 

1,12, 19H THRU ENERGY GROUP ,I2,5H )***) 

103 FORMAT! 1HK, 18H DELAYED SPECTRUM , I 2 , / ( 7 El 8 . 7 ) ) 

104 FORMAT! 71 10) 

C 

I F( KBEFF.EU.O) GO TO 1000 
NIT=NI J*NG 
DO 900 I J K =1 , N I T 
XN A ( I J K ) = XNA ( I J K) /XKEFF 
900 CONTINUE 

DO 1006 K J = 1 , N I J 

C THIS PORTION OF YJA IS USED TO STORE REGULAR PRODUCTION FOR 

C BETA EFFECTIVE CALCULATION 

1006 Y JA ( K J )=F { KJ ) 

GO TO 140 0 
1000 DQ 1005 I G= 1 , NG 
P( IG)-0.0 
S( I G) = C. 0 
A ( IG)=C.O 
R ( I G ) =C. 0 
SIXUG ) = 0.0 
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TL 1 ( I G ) = 0 * 0 
TL 2 ( IG ) = 0. 0 
Dlt IG ) =0 • 0 
1005 021 IG ) =0 • 6 
NG WP 1 = C 


THE FISSION PERTURBATION SOURCE 


IF (KP) 1 100,1100,1010 
1010 00 1025 I J= i , N I J 
1025 Ft I J ) = C • 0 

DO 10 5 C I G = 1 , NG 
K 1= IND EXt NTYP S, IG ,2) 

I J G = I NOE X ( N I J , I G , 0 ) 

DO 10 5 C I J = 1 , N I J 
IJG=I JG+1 
K=MA( I J) 

IF ( K ) 1050, 1050,1 030 
1030 ITEMP = INDEX(NTYNG , K , K1 ) 

Ft IJ )=F( IJ )+Ct I TEMP)*XN( I JGI 
1050 CONTINUE 


DO 1075 IG=l f NG 
IJG=INDEX ( N I J,IG,0) 

DO 10 7 C I J = 1 , N I J 
1 J G = I J G+ 1 

107 0 PI IG) =P( I G ) + XNA < I JG)*F < I J)*V< 1 J) 
107 5 PI I G ) =-CH I ( IG)*P( IG)/(PAP*XKEFF) 


THE SCATTERING PERTURBATION SOURCE 


N0wPl=-2 
NQWP 1 =— 1 
NUWP 1= 0 
NO wP 1 = 1 
N0wPl=2 


P-1 CROSS SECTIUNS**NZ.GT.l 
P-1 CROSS SECTI UNS**NZ. LE. 1 
DIFFUSIUN APPRUXI MATI UN 

P-0 TRANSPURT CORRECTED CRUSS SECT IONS **NZ .LE . 1 
P-0 TRANSPURT CORRECTED CROSS SECT I ONS **NZ .GT . 1 


1100 IF IKS) 1350,1350,1105 
1105 IF I NG— 1 ) 1350,1350,1110 

1110 DO 130 C I G = 1 , NG 
RSET=1 

DO 1115 I J= 1 , N I J 
1115 F ( IJ)=C.O 


UP SCATTERING 

IF I N UP I 1200,1200,1125 
1125 IF ( I G-NF AST) 1200,1130,1130 
1130 IF (IG-NG) 1135,1200,1200 
1135 K SE T= 2 
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IGG^IG 
I TY =N UP + 5 
1140 I GG-I GG+1 
I TY = I TY-1 

K 1= INDEX ( NTYPStIGf I TY) 

IJG=1NQEX (NIJ fIGGfO) 

C IJG IS THE GROUP SCATTERED FROM 

XJT=INDEX(NI J,IG,G) 

DO ii 5 C I J=i , NI J 
IJG=I JG*i 
K=MA( I J) 

IF (K) 1150,1150,1145 
1145 ITEMP = INDEXC NT YNG ,K ,Ki ) 

F( IJ)=F( l J)+C C ITEMP)*XNU JG) 

IF(NU WP1.GE.0) GO TO LL50 
Ijr = IJ T+l 

J MP U=N TYNG*NMAT+l TEMP 

IF( NO WP1. EQ.- 1) F( I J) =F ( I J ) +C ( JMPU) *XJ([ JG) 
IF( NGWPl.E0.-2) F( I J ) =F ( I J ) +C ( J MPU) * 

1< XJ( I JG)*XJA( I JT) + YJCI JG) *YJAU JT) ) 

1150 CONTINUE 

IF (IGG-NG) 1155,1200,1200 
1155 IF (ITY-5) 1200,1200,1140 

DU WN SCATTERING 

1200 CONTINUE 

IF C IG .LE. 1) GO TO 1275 
1205 K SET= 2 
IGG=0 

ITY = (NUP+4) + lG 
1210 1 GG= I GG+1 
I TY = I TY-1 

IF (NTYPS-I TY ) 1210,1215,1215 

1215 K1=INDEX(NTYPS,IG,ITY) 

IJG=INDEX(NIJ , IGG ,0) 

IJT=INDEX (Ni J, I G , 0) 

DO 12 5 C I J=1 , N I J 
I JG = I J G+l 
I J T = I J T+i 
K=M A ( I J) 

IF (K) 1250,1250,1245 

1245 IT£MP=IND£X(NTYNG,K,K1) 

IF( N0WP1.LT.0) GO TO 1246 

F{ IJ)=F(I J)+C (ITEMP)*XN(I JG) 

GO TO 12 5 0 

1246 JMPD=N rYNG*NMAT+I TEMP 

IF(NOWPl.EU.-l) F( IJ) =F { I J ) +C (JMPD) *XJU JG) 

I F { NO WPi.EQ.- 2) F ( I J ) =F ( l J ) +C < JMPD ) * 

I(XJ( IJG)*XJA( IJT)+YJ(I JG)*YJA(I JT) ) 

1250 CONTINUE 

IF< IGG-C IG-1) ) 1210,1280,1280 
1275 IF(KSET.EU.l) GO TO 1300 
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1280 IJG=INDEX(NIJ , IG , 0) 

C IJG IS THE GROUP SCATTERED TO 

IFINOWPl+l) 1283, 1282 , 1281 

1281 00 1285 I J=1,NI J 
UG = I JG+1 

S( IG) = S( I G) + XNA (I JG)*F(I J)*V( IJ) 

1285 CONTINUE 
GO TO 1301 

1282 DO 1286 I J=1,NI J 
IJG— I JG+1 

1286 S1X( I G )=S1X( I G ) +X JA ( I JG ) *F U J)*VU J) 

GO TO 1301 

1283 DO 128? I J = 1 * N I J 

1287 S1X( IG )=S 1X( IGH-F ( I J)*V(I J) 

C 

1301 IFIKAPROX.LE. l.OR.NOWPiiLT.O) GO TO 1299 
I F ( IG.NE.NG) GO TO 1300 
NQWP 1=— 1 

IFINZ.GT. 1) NOWPl=-2 
GO TO 1105 

1299 S( IG ) =— ( S { IG ) - SIX ( I G ) ) /PA P 

C TDSN USES P-1 CS WITH A FACTOR OF 3 INCLUDED 

1300 CONTINUE 


ABSORPTION 

*** 


SC ATTER IN G(REMOVAL) 

*** 

PERTURBAT ION 

LEAKAGE ( TRANSVERSE ) 

*** 


LEAKAG E( TRANSPORT) 

* * * 

LOSSES 


WITHIN GROUP SCATTERING*** 


1350 IF (KARl 1400, 1400, 1355 
1355 DO 138C IG=1,NG 

K 1= INDEX( NTYPS, IG ,1) 

I J G = I INDEX ( N I J , IG, 0) 

DO 13 75 I J = 1 , N I J 
IJG=I JG+1 
K=M A ( I J ) 

IF IK) 1375,1375,1370 

1370 I TEMP = INDE X(N TYNG »K,K1 ) 

L T E M = N TYNG*NM A T+NG* ( K- 1 ) + I G 

IF< KAPROX .EQ.2 ) L TEM=L Tt M+ N T YNG*NM AT 

LLTEM=LTEM+NG*NMAT 

TEMP=XNA< IJG)*XN( I JG)*V( I J) 

I F ( K.A P ROX . EQ . 0 ) GO TO 1371 
TEMPJ = 3.*XJA( I JG) * X J( I JG) * V ( I J) 

I F ( NZ .LE. 1) GO TO 1374 

TEMPJ = TEMPJ + 3.*YJA(I JG)*YJII JG)*V(I J) 
1374 CONTINUE 

IFIKAP ROX. EQ.2) TEMP=TEMP— TEMPJ 

1371 A( IG) =A( I G )+TEMP*C ( ITEMP) 

R ( IG) =R( I G)+TEMP*C ( I TEMP+3) 

TEMB=T£MP 

IF( KAPROX.EQ. 1) TEMB = TEMP- TE MP J 
TL II IG ) = TL1IIG) + TEMB*C(LTEM) 
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TL2( IG)=TL2(IG)+TEMB*C(LLTEM) 

S IGGG=CI I TEMP+2 )-CI I TEMP) -C < I TEMP+3 ) 
JTEMP=ITEMP+NTYPS*NG*NMAT 
SIGGG=SIGGG- { C € JTEMP+31 /3.0) 

IFIKAPROX.EQ. 1) 0 1 ( I G ) =D 1 « I G ) — TE MP J*C ( I T EMP+2 ) 

IF ( KA PRUX »E0 • 2) D1 ( IG ) =0 1 ( IG ) -TE MP J*S I GGG 

D1 IS THE TRANSPORT CROSS SECTION FOR KAPROX.EQ.l 

01 IS THE J-WGT G TOG CROSS SECTION FOR KAPR0X.EQ.2 

02 IS THE LEAKAGE CROSS SECTION FOR KAPROX . EQ.2 
IFIKAPROX.LE. 1) GO TO 1375 

D2( IG )=D2< IG)-TEMPJ*(C (I TEMP+2) -(C ( JTEMP+3 ) / 3 . 0 ) ) 
D2( I G ) = 0 2 I IG)+SIX( IG) 

1375 CONTINUE 

OH IG ) = D 1 ( I G ) /PAP 
U2( IG ) = D 2( IG ) /PAP 
TL1UG )=TL1(IG)/PAP 
TL 2 1 IG ) = TL2< IG) /PAP 
A ( IG) =A( I G )/P AP 
1380 K ( i G ) =R ( I G )/PAP 


1400 IF( KB EFF.EQ.O) GO TO 1800 

NOTE THE OELAYEO NEUTRON INFORMATION IS READ INTO OTHER ARRAYS 
< S AND YJA ) TO SAVE STORAGE 
READ! 5,101) ( S(I) , I=1,NDSP) 

SI I ) CONTAINS THE DELAYED NEUTRON FRACTION 
I SP AN = IGL 0- IGH I + 1 
DO 17 1 C IJK=1,NDSP 
LST=IN0EX( ISPAN.I JK,1) + NIJ 
LSP=L SI+I SPAN-1 

THIS PORTION OF YJA CONTAINS THE DELAYED SPECTRA 
1710 READ! 5,101) ( Y JA I K) , K=LST , LSP ) 

WRITE! 6,102) NDSP, IGHI ,IGLO 
DO 17 2 C I J K= 1 , NL) SP 
LST = INDEX( ISP AN, I JK » 1 ) + NIJ 
L SP =L SHI SPAN-1 

1720 WRITE ( 6,103) I JK, ( YJA ( K) , K=LST, LSP) 

CALCULATE THE IMPORTANCE OF THE DELAYED SPECTRA 

NST=N I J+NDSP*I SPAN 
DO 17 5 C I D = l , ND SP 
NST=NS7 + 1 
Y JA ( N ST ) = 0. 0 
C THIS PORTION OF YJA CONTAINS BETA 

DO 173 C IJ=i,NIJ 
TEEM = 0 .0 

IGA = INDEX I ISP AN, ID, 1) 

IGA=I GA+N I J 

DO 1740 I G=IGHI ,IGLO 

LOCXN=INDEXINI J,IG, I J > 

TEEM=TEEM +YJA ( IGA )*XNA ILOCXN) *S(ID) 
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1740 I GA = I GA + l 

TEEM=TEEM*YJA < 1JJ*VU J) 
1730 YJACNST )=YJA(NST)+TEEM 
YJA(NST) = YJA( NST) /PAP 
1750 CONTINUE 
C 

1800 RETURN 
C 

END 


$ I6FTC PER8 I 

SUBROUTINE DIF USE ( MA , V »DE LR , DE LZ , C , XN , X NA , 01 , D2 , T L 1 , T L2, A ) 


COMMON 

/C ALL1 / 


NG, 

NTYPS, 

1 NTYNG, 

NFAST, 

NUP, NR t 

NZ, 

NI J, 

2 

KF LUX , 

KBEFF « KCTYP, 

Ki N, 


3 

KP, 

KS, KAR, 

KD , 


4 KAPR3X,NMAT,NMIl> f NL)SP 

« I 3HI 1 1 GLO » XI R » X l Z 




COMMON /CALL2/ PAP 

DIMENSION MA( 1) , VUI «DELK (1 ) «DELZm tC Cl) * XN<1> «XNA( 1 ) » Dl( 1 ) * D2( 1) 
DIMENSION TLlll) ,TL2(1) ,A<1) 

I ND EX ( LENGTH, INDEXO, l NDEXU -LENGTH* ( INOEXO-L ) + I NOE XL 


THE DIFFUSION PERTURBATION LOSS 


104 FORMAT ( 71 10) 

105 FORMAT! 1HK, 69H BOUNDARY CONOITiCNS FOR RIGHT, LEFT, TOP, AND BOTT 

10M BOUNDARIES ** ,412) 

1405 DO 1595 I G = 1 , NG 

K1=INDEX( NTYPS,IG, 3) 

I F ( IG .NE • l.QR • KU. NE.-l ) GO TO 1409 
READ! 5,104) KRBC , KLBC , K TBC , KB BC 
WR I TE ( 6 , 1 0 5) KRBC , KLBC , KT0C , KBBC 

FOR THE FIRST DIRECTION (MUST BE AT LEAST THREE MESH INTERVALS) 
THREE POINT VANDERMONDE SULUTION FOR FLUX DERIVATIVE 
THE DIFFERENCE IN GRADIENT ACROSS PARALLEL FACES OF A MESH 
INTERVAL IS USED AS THE NET GRADIENT OF THE INTERVAL 

1409 X I =X I R 

DO 149 C i = 1 , NR 

I F ( I «E Q« 1 ) 0X1=2. *XI 

IF ( I-i) 1410,1410,1415 

1410 11=1+1 

GO TO 142 0 
1415 11=1-1 

DXI=2.*DELR( I 1)-DXI 
IF (I-NR) 1420,1425,1425 
1420 DEL l-DELRCl)+DEUMm 
DEL 2=0 ELR ( I )*DELR<11) 

DEL 3= DEL i*DE L 2 
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1425 DO 149C J =1»NZ 

I J = INDEX! NR, J , I ) 

K=MA( I J) 

IF (K) 1490, 149C, 1430 

1430 ITEMP = IN0EX(NTYNG»K»K1) 

IJG=INUEX(NIJ,IG,I J) 

IF (1-1) 1435,1435,1440 

1435 I JG4=I JG 

I JG5=I JG+ 1 
I JG6=I JG+ 2 
GO TO 1475 

1440 IF ( I— NR ) 1445,1450,1450 

1445 I JG4=I JG- 1 
I JG5=I JG 
IJG6=I JG+1 
GO TO 1475 
1450 I JG4=I JG- 2 
IJG5=I JG- 1 
IJG6=I JG 
C 

1475 12=1 

IF( I.EQ.NR) 12=1-1 

IF( l.NE.l.ANU.I.NE.NR) GO TO 1485 

I F ( KD . GE . 0 ) GO TO 1485 

IF( I2.NE. 1.3R.KLBC.EQ.0) GO TO 1479 

S IGN= 1.0 

NON E= I JG4 

N TWQ= I JG5 

GO TO 1480 

1479 IF C KRSC.EQ.O) GO TO 1480 
S IGN=- 1.0 

NUNE= I JG5 
N T WQ= I JG6 

1480 DEN0M=SIGN*DELR(12)*(DXI+DELR<I2)) 

A 2= ( XN l NT WO ) - XN ( NONE ) ) /OENOM 
A2A=( XNA( NTWO )-XNA(NONE) ) /OENOM 
0ELP=2.*SIGN*A2*0XI 
DELPA=2.*SIGN*A2A*DXl 

GO TO 1486 

1485 I F ( i.EQ.l) 12=1+1 

DEL P= 2 .* ( DELR { I 2 1 *XN ( I JG4) -DE LI * XN I I JG5 ) +D£LR( 12-1 )*XN( IJG6I )*OXI/ 
1DEL 3 

DEL PA =2.* I DELR ( 12 )*XNA ( 1 JG4 ) — 0£ L 1*XNA ( 1 JG5 ) + DELR( 12-1 l*XNA( IJG6) >* 
1DXI/0EL3 

1486 Dl( IG ) =0 1 ( 1G)-DELPA*C(I TEMPI* DELP*V(IJ) 

C SUCKLING LOSS - FLUX GRAOIENT CONTRIBUTION 

NGM=NG*NMAT 

LTEM=N TYNG*NMAT+NG*(K-1)+IG 
LLTEM = LTEM+NGM 
LLL=LTEM+ 2*NGM 

LLL IS LOCATION OF F I C K, S LAW CURRENT FACTOR 1/ ( 3*( SIGTR**2> ) 

LTEM AND LLTEM ARE LOCATIONS OF THE BUCKLING LOSS CROSS SECTIONS 
TL1( IG )=TL1< IGJ-C (LTEM) *D£ LP*OE LPA*C ( LLL) *V( I JI/PAP 
TL2( IG ) = TL2( IG)-C ( LLTEM) *DELP*DE LPA*C < LLL 1 *V ( I J )/PAP 
IFIKAR.NE.2) GO TO 1490 
I TM 1=1 TEMP-2 

A ( I G ) =A( I G )-C ( ITM1)*UELP*DELPA*C(LLL) *V< I JI/PAP 
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o n o o o n n 


C 

1490 CONTINUE 

Dl( I G ) =D 1 ( IG ) /PAP 
SEE EQ • 32 IN HRI TEUP 

FOR THE SECOND DIRECTION (MUST BE EITHER ONE OR THREE OR MORE MESH 
INTERVALS) THREE POINT VANDERMONDE SOLUTION FOR FLUX DERIVATIVE 
THE DIFFERENCE IN GRADIENT ACROSS PARALLEL FACES OF A MESH 
INTERVAL IS USED AS THE NET GRADIENT OF THE INTERVAL 

IF ( NZ . LE * 1) GO TO 1595 
XII=XIZ 
1500 DO 15 9 C J =1? NZ 

IF(J*EQ«1J DXII=2.*XII 
IF (J-l) 1505,1505,1510 

1505 J 1=J+1 

GO TO 1515 
1510 J1=J-1 

DX1 I=2**0ELZ( J 1 )-D XI I 
IF (J-NZ) 1515,1525,1525 
1515 DEL 1=DELZ(J)+DELZ(J1) 

DEL2=DELZ(J)*DELZ(J1) 

DEL 3=D£L 1*DEL 2 
C 

1525 IJ=INDEX(NRt J,0) 

IJG=INDEX(N1J,IG, I J) 

DO 15 9 C I =1, NR 
I J = I J + 1 
I J G= I J G+ l 
K=M A ( I J) 

IF (K) 1590,1590,1530 
1530 I TEMP = IND EX(N TYNG »K, K1 ) 

IF (J-l) 1535,1535,1540 

1535 I JG4= I JG 

I JG5= I JG + NR 
IJG6=I JG5+NR 
GO TO 1575 

1540 IF (J-NZ) 1545,1550,1550 
1545 1 JG 4= I JG- NR 
I JG 5= I JG 
I JG 6= I JG + NR 
GO TO 157 5 
1550 I J G 5= I JG-NR 
I JG4= I JG5-NR 
I J G6= I JG 
C 

1575 J 2 = J 

IF(J«EQ*NZ) J 2 = J- 1 

IF( J.NE.1*AND. J.NE.NZ) GO TO 1585 

IF(KD.GE.O) GU TO 1585 

IF{ J2.NE. l.OR.KBBC.EQ. 0 I GO TO 1579 

S I GN= 1*0 

NON E= I JG4 

NTWO= I JG5 

GO TO 1580 

1579 I F ( KT BC. E Q • 0 ) GO TO 1580 
S I GN — — UO 
NON E = I JG5 
N TWQ= I JG6 
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15 80 UENOM= SIGN*DEL2{ J2 )* <0X1 l +DE LZC J2) I 
A 2= ( XN ( NT WO ) - XN < NONE ) ) /DE NOM 
A2A = ( XNA( NTwQ )— XNA ( NONE ) ) /DE NOM 
DEL P = 2 .*S I6N* A2*D XI I 
DELPA=2.*SIGN*A2A*DX£ I 
GO TO 158 6 

1585 IHJ.EQ.li J2=J+i 

DEL P= 2 .* ( DELZC J2)*XNU JG4 l-DE LI * XN C I JG5 I+DELZ ( J2-1 ) *XN ( IJG6) )*DXII 
1/DE L 3 

DEL PA = 2** <DELZ( J2)*XNACI JG4 > -DE L 1 *XNA < I JG5H-DELZC J2-1 ) *XNA( IJG6H* 
1DXI I/DEL 3 

1586 D2( I G ) =D2 ( £ G ) -DELP A*C C £ TEMP) *DE LP* V ( I J) 

L TEM=N TYNG*NMAT+NG*(K-1) + IG 

LLL =L T£M+2*NGM 

TLK1G )=TL1( IG J-C < LTEM) *OELP*DELPA*C C LLL) *V < l J ) / PAP 
I F { KAR ,NE *2 ) GO TO 1590 
I TM 1= I TEMP-2 

A ( I G ) = AC I G )-C { I TMl)*DELP*DELPA*C(LLL)*V(IJ)/PAP 

C 

1590 CONTINUE 

D2C IG )=D2( IG) /PAP 
C 

c 

1595 CONTINUE 
RETURN 
END 


$ I8F TC PER9 

SUBROUTINE OUTPUT ( P, S, A, R, 01, D2 , Y J A, X KEFF , T LI , T L 2 ) 
C 

c 

c 

c 

c 

C THE COMMON STATEMENTS 


COMMON 

/CALL 1/ 


NG, 

NTYPS , 

1 NTYNG, 

NFAST* 

NUP, NR, 

NZ, 

NIJ, 

2 

KFLUX, 

KBEFF , KCTYP, 

KIN, 


3 

KP, 

KS, KAR, 

KD, 


A K APRQ X, N MAT, 

NMIO.NUSP 

, IGHI , I GLO ,XIR,XIZ 



COMMON 

/CALL2/ 

PAP 



COMMON 

/CALL3 / 

KXTKA1, KXTRA2, 

KXTRA3, 

KXTRA4, 

1 KXTRA5, 

EXTRA1 , 

£ XTRA2 , EXTRA3 , 

EXTRA4, 

EXTRA5 

THE DIMENSION 

STATEMENTS 



DIMENSION 

P (1* , 

S (11, 

AC 1 ) , 


1TL 1 C 1) ,TL2( 1J 

, R( 1> , 

D1 Cl I , 

D2C1) 

, YJAC 1) 


C THE FORMAT STATEMENTS 

100 FORMA T { 1H 1 ) 

10 L FORMAT (1HK,25H DELAYED SPECTRUM NUMBER ,I2,21H **** SPECTRAL SUM = 
1 ,F10.7/42X» 7HBETA = , F 1 0 .7 /37X , 12H BET A { EF F) = ,F10.7) 

113 FORMAT (7E16.7) 
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non 


114 F0RMAT(lHL*3t)H UNPER TURBED MULTI P LI C AT l ON FACTOR ****,F12.8/ 

139H PERTURBED MULTIPLICATION FACTOR **** , F 1 2 . 8 ) 

115 FORMAT! 1HK.4X ,10HTRANSVERSE ,6 X, 1 OHTR ANSVERSE , 56X ,6 HGROUP / 

15X, 10HL6AKAGE 1 , 6X, 10HLEAKAGE 2) 

116 FORMA T12E16.7 ,48X*E16.7) 

119 FORMA T ( 35X , 41H (( l /KEFF Pt RTUR BE 0 >- ( 1 / KEFF UNPERTURBED))) 

120 FORMAT ( 44X , 23HDI FF USI ON APPROXI MATI 0N//5X , 

11 1H PRODUCTION ,5 X,1 OH SCATTERING ,6X , 10 H A6S ORPT I ON, 6X ,7 HREMOV AL , 5X, 
216H DIF F. COEFF. 1 , 1 6H D1FF. COEFF. 2 ) 

121 FORMAT (1HJ.14H GROUP TOTALS) 

122 FORMAT (37X,37HTRANSP0RT CORRECTED P-0 APPROX I MAT ION// 5X , 

11 1H PRODUCTION ,5 X,1 OH SC AT TER I NG ,6X , 1 OH A 8S ORPT I ON, 6 X , 7HR EMO V AL , 9X , 
29H TRAN SPURT, 1 8H DIFF. COEFF. 2 ) 

123 FORMAT (42X,2 7HP— 1 TRANSPORT APPROX I M AT I 0N//5 X , 

11 IH PRU0UCTIJN,5X,10HSCATTERI NG ,6X , lOHABSORPT I 0N,6X .7HREM0VAL , 6X, 
214HG-T0-G***J-WGT,5X, 7H LEAKAGE) 


1000 WRITE (6, 100) 

IF! KBEFF.EQ.O) GO TO 1010 
l SP AN = IGLO- IGHI+l 
N ST =N I J +N D SP* I SP AN 
L U B =N I J 

DO 1015 ID = 1,NDSP 
NST=N ST + 1 
S UM =0 . C 

DO 102C I GD=1 , I SPAN 
LOB=LO B+l 

1020 SUM =SUM+YJA! LOB ) 

1015 WR 1 TE ! 6, 1 01 ) ID ,SUM, S! ID) , Y JA (NST) 

GO TO 110 0 
101C WRITE ( 6, 1 19) 

IF! KAPROX.EU. 0) WRITE(6,120) 

IF! KAPRUX.EQ. 1) WRITE(6,122) 

IF! KAPKOX .ED. 2) WRITE!fa,123) 

R HO =0 . 0 
R HOP =0.0 
R HQ S= 0 .0 
RHU A=0 .0 
RHQR = 0 .0 
RH0D1=C.0 
R HO 02 = C .0 
RHOTL 1 = 0. 0 
RHU TL 2 = 0. 0 
DO 1005 I G= 1 , NG 
RHO A=R HOA +A! IG) 

I F ( KAR .EQ .2) GO TO 1005 
RHUP = RHUP+P! IG) 

RHO S=RHOS + S( IG) 

R HUR=R HOR +R( IG) 

RH0D1 =RHO D 1+D 1 ! IG) 

RH002=RH0D2+D2! IG) 

1005 WRITE 16,113) P < IG ) , S l I G ) , A < I G) , Rt I G ) , D1 1 I G) , D2 ( I G ) 
WRITE (6,121) 

WRITE 16,113) R HOP, RHUS, RHO A , RH OR , RH0D1 , RH0D2 
WR ITE I 6,1 15) 

DO 10 5 C I G=1 , NG 

RHUG=A ( I G ) +R ( IG ) +P ( I G ) +S( IG)+TL1 (IG) +TL2! IG) 
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IF( KAPRUX.EQ. 0) RH0G=RHQG+D1 ( IGH-D2 ( IG) 
I F ( KA P ROX • GE . 1 ) RH0G=RH0S+01 ( IG) 

RHQ TL 1=RH0 TL i+TLl ( I G ) 

RHO TL 2=RHU TL2 + TL2 C I G ) 

RHO=RHO+RHOG 

1050 WRITfc<6,116) TL1( IG) , TL21IG) ,RHOG 
WRI TE( 6, 121) 

WR I TE { 6» 1 16) RH0TL1»RH0TL2»RHQ 
XKPER T=XKEFF/( 1.0+XKEFF*RH0) 

WRITE! 6» 114) XKEFF , XKPERT 
C 

c 

c 

1100 RETURN 
C 

ENU 


OVERLAY ORIGIN CARDS AND ASSIGNED LINK NUM8ERS 


$0R I G IN 

001 

IS 

LINK 

1* 

PARENT 

LINK 

IS 

0 

tOR I G IN 

□ 11 

IS 

LINK 

2 , 

PARENT 

LINK 

IS 

1 

SORiGIN 

OIL 

IS 

LINK 

3 , 

PARENT 

LINK 

IS 

1 

SDR i GIN 

Oil 

IS 

LINK 

4, 

PARENT 

LINK 

IS 

I 

SORIGIN 

Oil 

IS 

LINK 

5 , 

PARENT 

LINK 

IS 

1 

SORIGIN 

Oil 

IS 

LINK 

6, 

PARENT 

LINK 

IS 

1 

SORIGIN 

001 

I s 

LI NK 

7 f 

PARENT 

LINK 

IS 

0 

SORIGIN 

001 

IS 

LINK 

8 , 

PARENT 

LINK 

IS 

0 

SORIGIN 

001 

IS 

LINK 

9 , 

PARENT 

LINK 

IS 

0 
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* MEMORY MAP * 


SYSTEM 

FILE BLOCK ORIGIN 
FILES 1. UNIT06 

2. UNIT05 

PRE-EXECUTION initialization 
CALL ON OBJECT PROGRAM 
OBJECT PROGRAM 


OOOOO THRU 02717 
02720 

(NO BUFF POOL ATTACHED! 
(NO BUFF POOL ATTACHED! 
02750 
02757 

02764 THRU 16004 


L INK 

DECK 

ORIGIN 

CONTROL SECTIONS 

(/NAME/ 

= NON 0 LENGTH, ( LOC ) = DEL ET ED 

, *=NOT 

REFERENCED) 

0 

PER SN 

02764 

/// 

/ ( 1 7 1 30 ) 

/ CALL 1 / 

02765 

/CALL 2 / 

03015 

/CALL 3 

/ 03016 

/CINPT i 





03470 * 









.L INK 

03514 

/.LOT 

/ 03514 

/. LRECT/ 

03526 

/ .LV EC / 

03550 





-LXCON 

03572 

.LXSTR 

03572 

. LXSTP 

03575 

• L < OUT 

03617 * 

.LXERR 

03626 

. LXCAL 




.LXRTN 

03631 * 

IBEX IT 

0 36 31 

•DBCLS 

04011 * 

. L XAR G 

04100 * 

.LO 




.ELSE 

04115 * 

. LFBL 

04116 * 

• L J NB 

04117 * 

.OFUUT 

04120 



-LUVRY 

04124 

. LUVRY 

(04124) 

.LOT 

(03514 ) 

.LRECT 

( 0 35 26) 

. L VEC 

(03550) 



TGDM 

04635 

TAGDUM 

04635 









-LXSL 

04636 

. LXSEL 

04636 

. LXSL 1 

04637 

-LXTST 

04642 

.LXOVL 

04702 

• LXMOD 




-LXIND 

04770 

. LXDIS 

04773 

.LX FL G 

04774 

.L TCH 

04775 



.F p r rp 

05003 

-PEP T. 

05003 * 

. F Xfc M. 

05005 

.FX EM 

05005 * 

-STAG 

05234 * 

•FMCRT 




.CGUTU 

05363 * 

• FXUUT 

05367 * 

. FX AR G 

05401 * 

.EXIT. 

05403 * 

EXIT 




• D VCNK 

05410 * 

UVFLUW 

05411 * 

SYSONE 

05436 

.NOP 

05437 

.BLANK 


• ERAS. 

05467 

E. 1 

05467 

6.2 

05470 

E . 3 

05471 

F . 4 

05472 



• XLC . 

05473 

CC. 1 

05473 

CC.2 

05474 

CC. 3 

05475 

CC .4 

05476 



FCN V 

05477 

. FC UN • 

05477 

. DDPRE 

05507 * 

• FCNV . 

05515 

.ENOF S 

05526 

.CNVSW 




• F 0 X 1 

05534 

. F OX 2 

05535 

. DB C 

05537 

.DBC 14 

05602 * 

. ST0P4 




• DBC20 

05 621 

. DBCIO 

05635 

.DBC99 

05650 * 

. DUS W 

05673 * 

.FI XSW 




. D XPSt 

05716 * 

.DOZET 

06023 * 

. IC LO 

06064 * 

. STC1PJ 

06101 * 

.F C OUT 




. F CARD 

06 LI 7 * 

. ALCUD 

06121 * 

. AN P T 

06132 

.ONP T 

06147 

. L N T P 




.AOUT 

06262 

. UUUT 

06274 

.LOUT 

06323 

.GOUT 

06333 

. TG OUT 




.FL r 

06442 

.FXFL1 

06577 

.FX D 

06603 

.FXFL2 

06607 

.FXFL3 




. I NTS 

0661 7 

.TUPAC 

066 *5 

-W I UTH 

06641 

. FPACK 

C6646 

. T6 ST 




.DO SF 

06 7 00 * 

. KUUNT 

06701 

-LIST 

06704 

. DON E 

06713 

. mjTRF 




.3 MAR 

07 301 

.DEXP 

07 310 * 

.TEN 

07314 * 

. FBDBF 

07316 

.DATUM 




. WORD 

0 733 7 

. MUD 

0 7 340 

.Pt X 

07341 

.Ft XP 

07342 

.DIG 


F IOh 

07356 

.F 1 JH. 

0 7356 

. UDI U 

0 7403 * 

. DDL LT 

07544 * 

.UCPT 1 

0 7 747 * 

. F XP T 




.FF 1L. 

10176 

. DOF I N 

L0214 * 

• FR T N . 

10223 

. DDR TN 

LC227 * 



FIDS. 

10 3 32 

..F IUS 

10332 

.FCLS 

10354 * 

. .F IOC 

10401 

. . F Sf L 

10422 

. .fbck 




. . F rc K 

10443 

. IUEF 

10527 * 

..X EM 

105 35 * 

. . FC HK 

10552 * 



FROU . 

10665 

..FROU 

10655 









FrtRO . 

10 702 

. .F WHO 

10702 









FBuU . 

10730 

. . F BCD 

10730 

. .FBCW 

10732 

. .F BCB 

10743 





UN 1 T 06 

10 773 

. . UN 06 

10773 









UN 05 

10 774 

. UN 05. 

10774 









• lUE . 

10 7 75 











.IDE 66 

1100 7 











.R wUUfc 

11016 











.BCR EA 

11 035 

BCR t AO 

11035 

BCREAD 

( 11035 ) 







.BCR wD 

U 116 

. .BLRJ 

It LL5 

. . BC WO 

i 1 L17 * 

..BRD8 

11131 





FXP2 

IL167 

. XP2. 

11167 









.XtXP . 

11262 

. N u J E F 

11262 

. XEXP 

1 1266 * 

• Z J Ot F 

11271 





// 

17130 










1 

PER 1 

11274 

/// 

/( 1 7L30) 

/CALL! / (02765) 

/ C AL L 2 /( 030 15 ) 

/CALL 3 

/ ( 0301 6 ) 

/CINPT . 




/-HANS 

/ 1L275 

( NPUT 

1 3222 






2 

PER 2 

13246 

/CALL1 

/ ( 02 765) 

VRZ 

1 4 24 L 







PER 2 l 

14341 

RZREAO 

14500 








3 

PER 3 

13246 

/LALL1 

/(02 76 5) 

/CHANG /( 1 1275 ) 

CRE AD 

15673 





PER 4 

15 7o6 

XLHANh 

15776 








4 

PER 5 

13246 

I OAC A V 

L372U 








5 

PEK6 1 

L 3 246 

RE AD 1 T 

L 3450 









PERo 

13 516 

/CALLL 

/ (02 765 ) 

NREAD 

14330 







PER666 

14413 

/CALL1 

/( 02 76 5) 

REED 

14637 






6 

PER 7 

13246 

/CALL1 

/ ( 02 765) 

/CALL2 / (03015 ) 

PAP CAL 

13624 




7 

PER 8 

11274 

/CAL LI 

/(02 76 5) 

/CALL2 / (03015 ) 

/CALL3 / ( 0 30 16 ) 

PERTUR 

13454 


8 

PER 8 1 

11274 

/CALL1 

/ ( 02 765) 

/CALL 2 / (03015 ) 

OIFUSt 

13112 




9 

PER 9 

112 74 

/CALL1 

/ (02 765) 

/CALL2 / (03015 ) 

/ CAL L3 / ( 0 30 16 ) 

OUTPUT 

12200 


UNUSED CURE 



16005 

THRU L 7 1 2 7 







BEGIN EXECUTION- 


03030 


03S3I 
04107 * 


04744 

05350 * 
05403 * 
05452 


05530 
05605 * 
05702 
05 1 1 L * 
06223 
06342 * 
06413 
05 54 7 
0 704 7 
07324 * 
07343 
10145 * 

10441 * 


( 03030) 
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